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ABSTRACT 

This  report  describes  the  development  of  a  method  programmed 
for  a  digital  computer  to  evaluate  one-dimensional  detonation 
parameters  in  g^aseous  media.  The  calculations  are  based  on  the 
assumptions  that  the  mixture  is  comprised  of  ideal  gas  const!- 
tuents  and  that  the  wave  process  is  govetned  by  equilibrium  end 
conditions. 

Details  of  the  analytical  method  of  solution  are  presented 
Including  the  IBM  7090  computer  program,  described  in  Fortran 
language.  Sample  input  and  output  data  with  instructions  for 
preparing  the  input  data  are  given. 

Finally,  computations  are  performed  for  a  number  of  hydrogen- 
oxygen  mixtures  from  H2+O2  to  3H2+O2  at  several  initial  conditions 
covering  a  range  in  pressure  from  0.1  to  760  mm  Hg  and  in  temper¬ 
atures  from  -180  to  +200°F  and  a  comparison  is  made  with  calcula¬ 
tions  of  other  Investigators  for  a  stoichiometric  mixture  initially 
at  NTP. 
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NOMENCLATURE 

Dimensional  Variables. 

-  sound  velocity 

-  defined  by  eq.  6.  1.  17 

Cp-  -  molar  specific  heat  of  iXh  species 

Cp  -  mass  specific  heat  of  gaseous  mixture  - 

f»  -  defined  by  eq.  6.  1.  13 

-  molar  Gibbs  free  energy  at  Itll  species 

5*  -  molar  Gibbs  free  energy  at  tth  species  at  one  atmosphere 

-  molar  enthalpy  of  tth  species 

h  -  enthalpy  per  unit  mass  of  the  gaseous  mixture 
m;  -  molecular  weight  of  tlb  species 
m  -  molecular  weight  of  gaseous  mixture 

-  number  of  moles  of  Ltb.  species 

D  -  total  number  of  moles  of  gaseous  mixture 

-  pressure 

<lj  -  mole  fraction  of  component  in  hypothetical  system 
(R,  -  universal  gas  constant 
R  -  specific  gas  constant 

S  -  molar  entropy  per  unit  mass  of  the  gaseous  mixture 
T  -  absolute  temperature 
U  -  particle  velocity  relative  to  the  wave 
V  -  volume  per  unit  mass  of  the  gaseous  mixture 
X^*  -  mole  fraction  of  ith  species 
Yi  -  chemical  symbol  denoting  the  ctfa  species 

Zj  -  iteration  parameters  defined  by  eq,  6.  1.  14 
Q  ,  -  stoichiometric  coefficients  for  iih  chemical  reaction 

-  enthalpy  change  for  the  i%h  chemical  reaction  defined  by  eq.  6.  2.  3 

-  Gibbs  free  energy  change  for  the  Lib  chemical  reaction  defined  by 

eq.  6.  1.  15 
"  specific  heat  ratio 

-  chemical  potential  of  tib  species 
SubscriptF 

1  -  denotes  initial  properties  upstream  of  the  wave 

2  -  denotes  final  properties  downstream  of  the  wave 
f  -  denotes  properties  based  on  frozen  composition 

e  -  denotes  properties  based  on  equilibrium  composition 
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Dimensionless  Variables: 


C  - 

•  C,/R, 

y-  - 

h/R.T 

^rn  - 

m^/m 

M,  - 

u,/a, 

>^12  ■ 

P  - 

•  -Pj/f. 

1^^  -  pressure  ratio  at  CJ  state 

-  pressure  ratio  at  von  Neumann  spike 

-  pressure  ratio  in  Hugoniot  curve  where  V  =  1.  0 

0  -Vt, 

V  =  Vj/Vi  -  Ue/u, 


V.  - 

v«- 

V„- 

Vs- 

A0- 

AP- 

CJ 

At- 

AP„.- 


V  evaluated  from  Hugoniot  equation  at  given  p  and  approximate 

V  evaluated  on  Rayleigh  line  at  given  P 

V  evaluated  on  Rankine- Hugoniot  curve  at  given  P 

V  evaluated  from  equation  of  state  at  given  P  and  © 
correction  to  6  at  specified  P  on  Hugoniot  curve 
correction  to  F?q- 

correction  to 
correction  to 
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1.  INTRODUCTION 

The  calculations  of  detonation  parameters  have  been  made  by  several 
investigators,  notably  Lewis  and  Friauf  (1),  Berets  et  al  (2),  Edse  (3), 

Wolfson  and  Dunn  (4),  Luker  et  al  (5),  Eisen  et  al(6),  Bollinger  and  Edse  (7), 
Barrere's  (8),  Gordon  and  Zeleznik  (9).  and  Bird  et  al  (10). 

The  distinguishing  features  of  the  various  analyses  are  summarized  in 
Table  1.  All  have  in  common  two  basic  assumptions,  namely  that  complete 
thermodynamic  and  chemical  equilibrium  is  established  immediately  down¬ 
stream  of  the  wave,  so  that  dissipative  effects  and  variations  in  chemical 
reaction  rates  can  be  disregarded,  and  that  the  constituents  of  both  reactants 
and  products  behave  as  perfect  gases. 

The  work  in  references  (1),  (2),  (4),  and  (5)  followed  similar  approaches 
and  in  particular,  specified  the  GJ  state  as  that  for  which  the  final  Mach 
number,  based  on  the  frozen  sound  speed,  was  equal  to  unity.  Lewis  and 
Friauf  made  use  of  this  condition  in  their  equations  and  proceeded  to  calcu¬ 
late  detonation  properties  by  determining  successive  corrections  to  the 
assumed  initial  values  of  the  independent  variables,  TJ  and  using  an  iterative 
technique.  Berets  et  al  (2)  used  the  same  method,  while  Wolfson  and  Dunn (3) 
and  Luker  et  al  (4)  performed  the  calculations  in  a  similar  manner  but  with 
and  p^as  the  independent  variables. 

It  was  subsequently  shown,  however,  that  the  correct  specification  of 
the  CJ  state  must  be  based  on  the  equilibrium  sound  speed.  This  problem 
was  circumvented  by  Edse  (3)  and  Bollinger  and  Edse  (7)  who  determined  the 
CJ  state  from  the  condition  that  the  corresponding  wave  Mach  number  is  a 
minimum.  An  initial  estimate  was  made  for  the  CJ  value  of  the  final  tem¬ 
perature  in  both  cases,  the  corresponding  properties  on  the  Hugoniot  curve 
evaluated^and  then  corrections  to  the  CJ  temperature  found.  This  is  essentially 
the  approach  used  by  Eisen  etal(6)  who  chose,  however,  as  the  independent 
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parameter.  The  latter  method  gave  double-valued  solutions  for  each  value 
of  Mj,  and  the  CJ  state  was  determined  when  both  solutions  converged, 
yielding  the  minimum  v;ave  Mach  number.  In  addition,  Eisen  et  al  (6)  were 
the  first  to  evaluate  the  final  equilibrium  Mach  number  as  a  check  on  their 
CJ  calculations. 

Barrere  (8)  developed  a  graphical  method  for  calculating  detonation 
properties  which  gave  results  in  good  agreement  with  those  of  Eisen  et  al(6). 
The  CJ  state  in  this  instance  was  defined  as  the  point  on  the  Hugoniot  curve 
at  which  the  entropy  was  a  minimum.  Incorporating  this  into  the  conservation 
equations,  1^  and  p^were  used  as  the  independent  variables  to  evaluate  the  CJ 
properties. 

Zeleznik  et  al  (9)  were  the  first  to  actually  utilize  in  their  computations 
the  principle  that  at  the  CJ  state  the  equilibrium  Mach  number  is  unity.  On 
this  basis  initial  estimates  were  made  for  the  independent  variables  and 
72  and  corrections  found  until  the  Hugoniot  equation  and  equation  of  state 
were  satisfied. 

Bird  et  al(lO)  computed  the  CJ  properties  by  finding  the  point  of  minimum 
wave  Mach  number.  Since  the  thermodynamic  properties  are  functions  of  the 
final  temperature,  ,  alone,  this  parameter  was  chosen  as  the  independent 
variable,  eliminating  the  need  for  re-evaluation  of  these  properties  at  each 
iteration.  The  program  is  generalized  to  include  the  presence  of  condensed 
phases  in  the  products  of  reaction. 

In  the  present  work,  a  method  was  developed  for  the  determination  of 
detonation  and  deflagration  parameters  including,  in  addition  to  the  CJ  state, 
solutions  on  both  the  deflagration  and  detonation  branches  of  the  Hugoniot 
curve.  The  procedure  was  programmed  for  an  IBM  7090  digital  computer  in 
Fortran  language.  For  detonations  the  calculations  are  terminated  when  states 
with  pressure  ratios  corresponding  to  the  von  Neumann  spike  on  the  one  end. 
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and  the  condition  that  the  density  ratio  is  equal  to  unity  on  the  other,  are 
attained.  For  deflagrations,  the  calculations  are  restricted,  of  course,  to 
values  of  the  pressure  ratio  equal  to  and  less  than  unity.  It  should  be  noted 
also  that  since  the  allowable  chemical  reactions  are  specified  in  the  input 
data,  the  parameters  for  other  processes,  such  as  non-reactive  or  dissoci¬ 
ative  waves,  may  be  easily  determined. 

During  the  initial  phases  of  this  study,  the  wave  velocity  was  chosen 
as  the  independent  variable,  following  the  approached  used  by  several  other 
investigators.  For  each  specification  of  the  velocity,  then,  solutions  on 
both  the  strong  and  weak  branches  of  the  tlugoniot  curve  were  obtained  and 
the  CJ  state  determined  when  the  solutions  converged  at  the  minimum  value 
of  wave  velocity.  However,  difficulties  in  obtaining  convergence  within  a 
desired  accuracy  were  encountered  and  it  has  been  then  found  that  the' most 
convenient  independent  variable  is  the  pressure  ratio. 

Since  the  thermodynamic  properties  of  each  constituent  are  uniquely 
functions  of  the  temperature  alone,  the  temperature  ratio  suggests  itself 
as  a  convenient  independent  variable.  However,  as  shown  in  Figs.  10  and  11, 
where  temperature  and  pressure  ratios  are  plotted  respectively  against  the 
dimensionless  entropy,  the  variation  in  the  pressure  ratio  is  about  ten  times 
greater  than  that  in  the  temperature  ratio  along  the  Hugoniot  curve,  so  that 
choosing  the  pressure  ratio  as  the  independent  variable  allows  greater 
latitude  in  making  an  initial  estimate  of  its  value  at  the  CJ  state. 

The  CJ  state  was  specified  by  the  condition  that  the  local  equilibrium 
Mach  number  is  equal  to  unity  which,  within  the  accuracy  of  calculations, 
corresponds  exactly  to  the  state  of  minimum  wave  velocity. 

Calculations  of  detonation  parameters  were  carried  out  for  several 
hydrogen- oxygen  mixtures  over  the  range  of  initial  conditions  as  follows: 
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Initial  Initial 


The  best  available  source  of  thermodynamic  data,  namely  the  JANAF 
Tables  (11)  prepared  by  the  Dow  Chemical  Company,  was  used  in  the  cal¬ 
culations.  Within  the  accuracy  of  these' data,  the  results  for  the  CJ  wave 
were  found  to  be  in  good  agreement  with  those  of  references  6,  8,  and  9 
obtained  for  similar  conditions.  In  addition  the  CJ  wave  velocity  exhibits 
an  increase,  independently  of  composition,  of  approximately  15  percent  as 
the  initial  pressure  was  increased  from  0.  1  to  760  mmHg,  and  a  decrease 
of  approximately  3  percent  as  the  initial  temperature  was  increased  from 
-180  to  +200^F.  Over  the  same  range  of  initial  conditions,  the  CJ  pres¬ 
sure  and  temperature  ratios  were  found  to  be  more  sensitive  functions  of 
tne  initial  temperature  than  of  the  initial  pressure. 
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2.  DETERMINATION  OF  THE  HUGONIOT  CURVE 

The  detonation  wave  is  considered  to  be  a  one -dimensional  front  and  the 
process  to  take  place  in  the  absence  of  dissipative  effects  so  that  the  simple 
continuity,  momentum,  and  energy  equations  apply.  The  conditions  of  the 
initial  gas  mixture,  including  its  composition,  are  assumed  to  be  known, 
while  the  thermodynamic  properties  of  each  constituent  of  both  the  initial 
mixture  and  final  products  must  be  specified  as  input  data.  Each  gas  species 
is  assumed  to  obey  the  ideal  gas  law  so  that  the  thermodynamic  data  of  each 
species  can  be  expressed  conveniently  in  the  form  of  a  polynomial  expansion 
in  terms  of  the  temperature.  Any  number  of  chemical  reactions  may  be 
considered  in  the  equilibrium  system,  including  the  case  of  no  chemical 
reaction,  corresponding  to  the  Rankine-Hugoniot  curve. 

The  final  conditions  of  the  wave  are  determined  from  a  solution  to  the 
conservation  equations  and  the  equilibrium  equation  of  state.  The  equations 
expressing  conservation  of  mass,  momentum,  and  energy  respectively  are: 


II 

2.  1 

2  2 

^  V,  Vj 

2.  2 

2.3 

where  u  is  the  particle  velocity,  V  the  specific  vobune,  the  pressure, 
and  h  the  enthalpy  while  subscripts  1  and  2  denote  initial  and  final  states 
of  the  wave,  respectively.  The  enthalpies  in  the  energy  equation  are 
evaluated  by  means  of  a  simple  summation  procedure: 


2.4 
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wiicre  the  ,  the  mole  fractions  of  the  species,  are  determined  for  state  1 
from  the  specified  composition  of  the  initial  mixture  and  for  state  2  from  the 
composition  corresponding  to  loo&l  thermodyna.'ic  equilibrium*  The  mean  mole- 
caiar  weight  is  a  weighted  ayerage  of  the  molecular  weights^  rh^  #  of  the 

,f 

species  given  by  the  relation! 

m  2. 5 

and  h-  ,  the  mclar  enthalpy  of  the  species,  can  be  expressed  as: 

» 

where  h.  is  the  enthalpy  of  formation  at  a  reference  temperature  "TJ  and 
is  the  specific  heat  of  the  species. 

The  conservation  equations,  eqs.  (2.  1),  (2.  2)  and  (2.  3),  can  be  com¬ 
bined  to  yield  the  well-known  Hugoniot  equation: 

^t)  ^ 

which,  expressed  in  dimensionless  terms,  has  the  form: 

-Mp-jXV+I)  2.S 

where  P  is  the  pressure  ratio  V  the  volume  ratio,  and  y  is 

the  dimensionless  enthalpy  defined  as: 

where  T  is  the  temperature  and  the  specific  gas  constant,  is  expressed 
in  terms  of  the  universal  gas  constant  ^  and  the  mean  molecular  weight: 


R=  (a/r 


2.  10 
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The  equation  of  state  of  the  madiuB  behind  the  mre  oan  be  written  as  t 

7^^^“  R2I2 

In  dimensionless  form,  eq.  (2.  11)  becomes 

i 

F\Jne/?7? 

where  ©  and  TH  are  respectively  the  temperature  and  molecular 
weighs  ratic3,  Vi  and  m^mi ,  across  the  wave. 

Since  it  is  most  convenient  to  work  in  the  specific  volume- 
temperature  plane,  eqs.  (2.8)  and(2.  12)  can  be  rearranged  to  yield 
the  following  expressions: 


v-zlVr^MF-i)  - 1 

2. 13 

V  *  &/(pm) 

2. 14 

Both 7^  and  TH  can  be  determined,  of  course,  once  P  and  Q  are 
known;  hence,  by  specifying  a  value  of  P,  solutions  to  sqs.  (2.  13) 
and  (2.  14)  can  be  found  as  shown  in  Fig  1.  The  intersection  of  the 
equation  of  state  (2.  14)  and  the  Hugoniot  equation  (2.  13),  which  repre¬ 
sents  a  solution  point,  i.  e.  ,  a  point  on  the  Hugoniot  curve  correspon¬ 
ding  to  final  conditions  of  the  wave,  is  found  in  the  following  manner. 

An  initial  choice  is  made  for  Q  and,  at  the  specified  P,  the  equili¬ 
brium  composition  and  corresponding  thermodynamic  properties  are 
caluclated.  If  the  initial  choice  of  Q  was  correct,  then  eqs.  (2.  13)  and 
(2.  14)  would  be  satisfied  simultaneously  within  the  prescribed  accuracy. 
If  not,  a  correction  to  G,  A0  ,  is  obtained  by  expansion  of  eqs.  (2.  13) 
and  (2.  14)  in  terms  of  a  Taylor  series  about  the  initial  approximation  to 
@.  Equating  the  resulting  expansions  and  retaining  only  first  order 
terms,  we  obtain  then; 


2.  11 


2.  12 


A0  =  -  (^)p,  J 


2.  15 
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where  andVj  are  the  volume  ratios  evaluated  from  the  Hugoniot  equation 
and  equation  of  state  lespectively  at  the  specified  P  and  previous  value  of  G. 
Differentiating  eqs.  (2.  1 3)  and .(2.  14)  with  respect  to  ©  at  constant  P  yields, 


respectively: 


AL. 

p-i 


i  lo^  3lp 


2.  16 


2.  17 


where  C-"^  is  the  non-dimensional  equilibrium  specific  heat  at  constant 
P  and  the  subscripts  H  and  S  denote  differentiation  in  the  V~©  plane  along 
the  path  of  the  Hugoniot  equation  and  along  the  path  of  the  equation  of 
state  respectively.  The  derivation  of  equation  (2.  17)  is  given  in  section 
6.  5.  Successive  corrections  are  made  to  the  temperature,  ©  ,  until  the 
values  of  V  computed  from  eqs.  (2.  13)  and  (2.  14)  agree  within  a  pre¬ 
scribed  accuracy, 

Eqs,  (2.  13)  and  (2,  14)  can  also  be  solved  by  choosing  0  as  the 
independent  variable  and  iterating  for  the  correct  value  of  P  , 

Eq,  (2.  15)  then  becomes 


A©  “  (Vff  “”V< 


2.  18 


where  the  subscripts  H  and  S  now  refer  to  differentiation  in  the  V~P  plane 
along  the  Hugoniot  equation  and  the  path  of  the  equation  of  state  respec¬ 
tively.  For  this  case  eqs,  (2.  16)  and  (2,  17)  become  respectively 


2.  19 


2.  20 
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£q.  (2.  19)  is  obtained  by  straightforward  differentiation  of  eq.  (2.  13), 
while  eq.  (2.  20)  is  derived  in  section  6.  5.  Fig.  2  shows  the  solution  for  a 
point  on  the  Hugoniot  curve  in  the  V’P  plane.  The  initial  conditions  are  the 
same  as  those  of  Fig.  1  and  the  value  of  @  ,  Fig.  2,  was  chosen  so  that  both 
soltrticm^  correspond  to  some  point  on  the  Hugoniot  curve.  A  comparison  of 
Figs.  1  and  2  shows  that  the  Hugoniot  equation  and  equation  of  state  are  nearly 
orthogonal  in  the  (constant  pressure)  plane  while  in  the  V’P  (constant 
temperature)  plane  they  are  almost  parallel.  Hence,  in  addition  to  a  wider 
latitude  for  the  initial  choice  of  the  independent  parameter,  choosing  P  , 
rather  than  S  ,  as  the  dependent  variable  should  minimize  problems 
associated  with  convergence. 

For  purposes  of  illustration.  Fig.  3  shows  the  Hugoniot  curve  for  a 
mixture  initially  at  60®F  and  1  atmosphere  in  the  P-V  plane.  Shown 
also  is  the  CJ  Rayleigh  line  whose  intersection  with  the  Rankine- Hugoniot 
curve  defines  the  von  Neumann  spike. 
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3.  DETERMINATION  OF  THE  CHAPMAN  JOUGUET  STATE 

*  • 

The  CJ  state  is  defined  as  that  state  at  which  the  final  equilibrium  Mach 
number  is  equal  to  unity.  It  is  determined  by  finding  the  intersection  of  the 
Hugoniot  curve  with  the  line  =  1 ,  in  the  H/P  plane.  Fig.  4  shows  the 

Hugoniot  curve  in  the  P  plane  for  a  mixture  initially  at  60^F  and 

1  atmosphere,  and  the  line  Me, a  *1.'  The  CJ  state  is  determined  by  making  an 
initial  guess  for  the  CJ  pressure,  ^7,  linearizing  the  Hugoniot  curve  in  the 
Mpp'P  plane,  and  solving  for  corrections  in  the  initial  choice  of  .  At 
each  approximate  value  of  Pc.-  ,  the  corresponding  properties  on  the  Hugoniot 
curve  are  evaluated  by  the  method  described  in  the  previous  section. 

While  previously  it  was  possible  to  find  analytical  expressions  for  the 
derivatives  of  the  governing  equations,  a  numerical  method  must  be  used  here. 
The  derivative  of  the  Hugoniot  curve  in  the'  P  plane  is  found  by  taking  a 

small  increment,  SP  ,  in  the  previous  estimate  of  F^j  and  evaluating  the 
corresponding  incremental  change  in  2  ^  SMe,?  .  Hence,  to  a  good  approxi¬ 

mation,  the  derivative  of  the  Hugoniot  curve  in  the  P  plane  is  given  by 
,  and  it  can  be  easily  verified  that  the  correction  -A  j  is  given  by 

APcr  =  3.1 

It  is  advisable  to  overestimate  the  initial  approximation  to  f?r  in 

order  to  insure  that  the  corresponding  value  of  \/  on  the  Hugoniot  curve  will 
be  less  than  one.  The  stability  of  the  iteration  procedure  is  improved  by  re¬ 
quiring  that 

1APct)<  1  3.2 

which  should  prevent  the  new  approximation  to  Pt  from  becoming  too  small 
and  yielding  an  imaginary  wave  velocity.  The  iteration  is  continued  until  the 
value  of  Me 2  is  sufficiently  close  to  one,  so  that  ARy  is  less  than  a  pre¬ 


scribed  value. 
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4.  DETERMINATION  OF  THE  VON  NEUMANN  SPIKE 


The  von  Neumann  spike  conditions  correspond  to  the  properties  behind 
a  non- reactive  shock  wave  which  propagates  at  a  Mach  number  equal  to  that 
of  the  CJ  detonation.  In  this  report  the  von  Neumann  spike  is  used  to  establish 
the  upper  limit  to  the  Hugoniot  curve.  The  state  in  the  P-Vplare  corresponding 
to  the  von  Neumann  spike  is  determined  by'  the  intersection  of  the  Rayleigh  line, 
which  emanates  from  the  initial  state  and  is  tangent  to  the  Hugoniot  curve,  with 
the  Rankine- Hugoniot  curve.  Fig.  3  illustrates  the  solution  for  the  case  of  a 

mixture  initially  at  60^F  and  1  atmosphere.  The  iterative  technique 
used  in  the  method  of  solution  is  similar  to  that  outlined  in  section  3.  An  initial 
guess  is  made  for  ,  the  von  Neumann  spike  pressure  ratio  and  the  equations 

of  the  Rayleigh  line  and  Rankine -Hugoniot  curve  are  linearized  permitting  cor¬ 
rections  to  be  found  to  the  initial  choice  of  Pvn  • 

It  is  necessary  to  know  the  value  of  the  derivatives,  ^  ,  along  both  the 

Rankine -Hugoniot  curve  and  the  Rayleigh  line  in  the  P’V  plane,  in  order  to 


linearize  the  equations  describing  these  curves.  The  first,  fp]  ,  is  found 
by  taking  a  small  increment  in  P  ,  S  P  ,  and  evaluating  the  corresponding 
increment  in  V  ,  SV  ;  the  ratio  ^  ,  yields,  then,  a  sufficiently  good 

approximation  to  the  derivative,  .  The  second,  invariant  and 

may  be  found  analytically.  The  dimensionless  Rayleight  line  equation  is  given 


,  ,  is  invariant  and 


(P-l)=  4.1 

where  ,  and  are,  respectively,  the  initial  specific  heat  ratio  and  Mach 
number.  Since  both  and  M,  remain  constant  along  the  Rayleigh  line, 
differentiation  of  eq.  (4.  1)  yields. 


4.2 
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An  initial  approximation  to  the  von  Neumann  :Lpike  pressure,  ,  is 
made  and  at  this  pressure,  the  required  derivatives  are  evaluated  as  outlined 
above.  It  can  easily  be  verified  that  the  correction  to  the  approximation  of 
PvN  is  given  by: 

where  and  \4  are,  respectively,  the  values  of  the  volume  ratio  evalu- 
ated  from  the  equations  for  the  Rankine-Hugoniot  curve  and  Rayleigh  line 
at  the  previous  value  of  The  iteration  is  continued  until  is  sufficiently 

close  to  so  that  is  less  than  a  prescribed  value. 
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5.  DETERMINATION  OF  THE  POINT  ON  THE  HUGONIOT  CURVE  AT  V  =  1.  0 

The  pressure  ratio,  Fv^/  «  corresponding  to  the  point  on  the  Hugoniot 
curve  where  V»1  is  used  to  establish  the  lower  limit  of  the  Hugoniot  curve. 

It  is  determined  by  the  intersection,  in  the^-P  plane,  of  the  Hugoniot  curve 
with  a  line  of  45^  slope  passing  through  the  ,  P*  0  ,  corresponding 

to  the  Hugoniot  equation  with  V-l.O,  The  two  equations  are  then: 


5. 1 

5.2 

Similar  techniques  are  employed  to  solve  eqs.  (5.  1(  and  (5.  2)  as  that 
used  in  sections  2,  3,  and  4.  An  initial  approximation  is  made  for 
the  two  equations  (5.  1)  and  (5.  2)  are  linearized  at  this  value  of  ,  and  a 
correction  found  for  the  initial  guess. 

The  derivative  of  eq.  (5.  1)  in  the  plane  is  found  by  taking  a 
small  inorement  S P  in  the  approximate  and  evaluating  the 

corresponding  increment  in  .  Hence  a  good  approximation  to  the 

derivative  of  eq.  (5.  1)  is  given  by  .  The  derivative  of  eq.  (5.  2)  is 

obviously  equal  to  one.  It  can  be  easily  verified  then  that  the  correction 
to  the  initial  approximation  of  is  given  by: 

5.3 

where  and  ^  F^re  evaluated  from  eqs.  (5.  1)  and  (5.  2)  respectively 
at  the  approximate  F^.,  . 

The  iteration  is  continued  again  until  is  sufficiently  close  to 

^  .  that  AP^.,  is  less  than  a  prescribed  value, 

H 
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6.  EVALUATION  OF  EQUILIBRIUM  COMPOSIHON  AND 
EQUILIBRIUM  THERMODYNAMIC  DERIVATIVES 

The  preceding  analysis  requires  the  evaluation  of  the  equilibrium  ccm" 
position  of  a  multicomponent  mixture  of  gases  at  a  specified  temperature 
and  pressure.  Brinkley  (12)  has  developed  a  convenient  method  for  equili¬ 
brium  calculations  for  use  in  machine  computation,  and  this  procedure  has 
been  adopted  for  our  work. 

6.  1  Equilibrium  Composition 

Following  Brinkley  (12)  and  Obert  (14)yWe  shall  consider  the  products 
to  consist  of  a  certain  number  of  species,  referred  to  as  constituents,  of 
which  a  certain  minimum  number  are  components.  The  distinction  between 
constituents  and  components  is  derived  from  the  fact  that  in  a  closed  svstem 
at  a  given  pressure  and  temperature  in  which  departures  from  equilibrium 
are  considered,  a  component  is  defined  as  an  independent  variable  while  a 
constituent  is  defined  as  a  dependent  variable.  (In  an  equilibrium  system 
at  a  given  pressure  and  temperature,  the  composition  is  fixed  and  herxce 
the  components  may  not  be  varied  independtly).  For  example,  consider  a 
closed  system  in  which  H2»  ,  and  H2O  are  present  at  a  given  pressure  and 

temperature.  Any  two  of  the  three  constituents  may  be  considered  as  com¬ 
ponents.  Formal  rules  for  selecting  components  were  given  by  Brinkley  (13). 

The  constituents  of  the  equilibrium  system  can  be  expressed  symbcli- 
cally  in  terms  of  the  components  in  the  following  form: 


where  designates  the  chemical  symbol  of  the  constituent,  Y  -  the 
chemical  symbol  of  the  jth  component,  and  are  the  stoichiometric  co¬ 
efficients.  In  this  section,  L  subscripts  refer  to  constituents  and  j  and  k 
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subscripts  refer  to  components.  For  example,  consider  the  hydrogen- oxygen 
system  in  which  six  constituents  are  assumed  present,  namely.  Hj  .  Os  .  HjO. 
OH  .  H  ,  and  O  ,  and  let  the  components  be  HsO  and  Os.  Then  the  matrix 

formed  by  the  s  is: 


— '-^..,Components 

*  Constituents 

HzO 

Oe 

Ht 

1.00 

-0.  50 

Oi 

0.00 

1.00 

H.0 

1.00 

0.  00 

OH 

0.  50 

0.  25 

H 

0.  50 

0.  25 

o 

0.  00 

0.  50 

The  composition  of  the  initial  mixture  determines  the  mole  fraction  of 
each  component  in  a  hypothetical  system  consisting  of  components  only. 
The  mass  balance  between  the  hypothetical  system  and  the  equilibrium  con¬ 
stituents  is  expressed  then  by  the  relationship: 

where  is  the  number  of  moles  of  the  constituent  in  the  equilibrium 
system  derived  from  the  hypothetical  system.  The  mass  action  relations 
for  the  system  under  consideration  can  be  expressed  as: 

where and/ij  are  the  chemical  potentials  of  the  constituent  and  jtli 
component  respectively. 
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Dividing  eq.  (6. 1.  2)  byn  >  where 


leads  to: 


6. 1.4 


6.1.5 


which,  when  summed  over  the  j  components,  becomes  after  rearrange¬ 
ment: 


6.1.6 


where  ^  is  defined  by  the  relation: 


Substituting  eq.  (6. 1.  6)  into  eq.  (6. 1.  5)  yields  finally: 


6. 1.7 


6. 1.8 


A  quantity  G-  can  be  defined  now  as  the  residue  of  eq.  (6.  1. 8),  i.  e.  , 

J 

where  Gj  is  identically  zero  when  the  correct  set  of  are  substituted 
into  eq.  (6.1.9). 

The  chemical  potential  of  the  tth  ideal  gas  constituent  may  be  written 
in  the  form: 

(fiTlo^ipx^  6.  1. 10 

where  is  the  chemical  potential  of  the  constituent  at  the  temperature 
in  question  and  one  atmosphere.  Note  that  the  chemical  potential  of  a  pure 
species  is  equal  to  the  Gibbs  free  energy  per  mole.  Hence,  for  this  case: 


17 


where  g?  is  the  molar  Gibbs  free  energy  of  the  constituent.  Substituting 
eq.  (6. 1. 10)  into  eq.  (6. 1.  3)  and  rearranging  terms  yields: 

where 

^  "  /og  (f)  ~ 

and  Z:  ,  the  iteration  parameter,  has  the  formal  definition: 

The  quantity.  A^g*  .  used  in  eq.  (6. 1. 13)  is  defined  by  the  relation: 

Now  an  initial  guess  must  be  made  for  .  Then  substituting  the 
value  of  X/from  eq.  (6. 1. 12)  into  eq.  (6. 1. 9),  expanding  the  resulting 
equation  in  a  Taylor  series  about  the  approximate  set  of  Sj  .  and  retaining 
only  first  order  terms,  one  obtains  a  set  of  linear  equations  of  the  form: 

where 

""  A  ^ 

and  Af,  are  the  correction  terms  for  the  initiaichoice  of  Zj .  With  the  new 
values  of  2  .  the  procedure  is  repeated  and  the  iteration  continued  to  the 


desired  accuracy. 
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/-Ai£#x4) 

U/qoTIp^^^  Ijloa-p/r 


6.2  Evaluation  of  thfe  Derivatives  y/f'^TA^  lj|nj-^/r 

Differentiating  eq.  (6.  1.12)  with  respect  to  logCD  at  constant  -f  one 

obtains:  _ _ 

log  Xt'  )  _  A:h  _ 


/ d  log  Xi  ]  -  A;h  ^  Q  j'^'Zj  \ 
wlogT  l^p  {RT  yP*'iwlo9'^' 


6.  2.  ] 


where  we  have  made  use  of  the  thermodynamic  relation 

-  A-h 

^  logT/^  T 


6.  2.2 


and  A^h  is  defined  as: 

6.2.3 

Taking  the  derivative  of  eq.  (6.  1.  8)  with  respect  to  logfr)  at  constant  'p , 
where  the  correct  set  of  have  been  used,  and  substituting  eq.  (6.  2.  1)  into 
the  resulting  expression  yields: 


6.  2.4 


which  gives  a  set  of  simultaneous  equations  from  which  the  set  of  derivatives 
may  be  determined,  allowing,  finally,  the  evaluation  of  the  set  of 

derivatives  from  eq.  (6.2.1) 

The  set  of  derivatives  derived  in  a  similar  manner. 

Differentiating  eq.  (6.  1.  12)  with  respect  to  \oq(^  at  constant  T  we  have  now: 


6.  2. 


while  taking  the  derivative  of  eq.  {6.  1.8)  with  respect  to  -p  at  constant 
T  ,  where  again  the  correct  set  of  Xi  have  been  used,  and  substituting 
eq.  (6.  2.  5)  into  the  resulting  expression  we  obtain  finally: 


6.  2.6 
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This  provides  a  set  of  simultaneous  equations  from  which  the  set  of  derivatives 
may  be  determined.  Consequently*  the  derivatives  (i 
may  be  evaluated  from  eq.  (6!  2,  5). 


6.  3  Partial  Derivatives  of  Molecular  Weight 

The  partial  derivatives  of  the  molecular  weight  can  be  expressed  most 
conveniently  in  terms  of  the  partial  derivatives  of  thr  mole  franc tions.  Making 
use  of  a  mathematical  identity,  we  can  write: 


.  J. 

'f  m 


6.3. 1 


With  the  use  of  eq.  (2.  5),  where  is,  of  course,  a  constant,  we  obtain; 

A1 


6.3.2 


Similarly,  it  can  be  shown  that: 


6.3.3 


6.  4  Evaluation  of  Equilibrium  Specific  Heat 

The  equilibrium  specific  heat  per  unit  mass  is  denoted  by  the  expres¬ 
sion  (^^ip  where 

h 


6.  4. 1 


Performing  the  indicated  differentiation  yeilds  then: 

where  the  first  term  on  the  right  hand  side  represents  the  frozen  composition 
specific  heat  per  unit  mass,  .  Expressing  the  derivatives  in  logarithmic 
form  we  have: 


6.4.3 
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Note  that  for  the  case  of  frozen  composition,  the  derivatives  on  the  right  hand 
side  of  eq.  (6.  4.  3)  vanish,  leading  to  the  identity  that  • 


6.  5  Evaluation  of  Thermodynamic  Derivatives  CaTa«  and  U-p, 


Differentiation  of  the  ideal  gas  equation  of  state  with  varying  molecular 


weight  yields: 


6.5. 


These  relations  are  necessary  to  evaluate  the  equilibrium  sound  velocity  and 
specific  heat  ratio,  while  eq.  (6.  5.  1)  has  already  been  used  in  section  2. 

The  derivatives  on  the  right  hand  side  of  the  equations  have  been  determined 
already  in  section  6.  3. 


6.  6  Evaluation  of  the  Equilibrium  Specific  Heat  Ratio 

The  equilibrium  specific  heat  ratio,  ,  is  defined  as: 


6.  6.  1 


which  upon  substituting  the  well-known  specific  heat  equation: 


6.  6.  2 


becomes: 


6.  6.  3 


Substituting  now  eqs.  (6.  5.  1)  and  (6.  5.  2)  into  (6.  6.  3)  and  rearranging  terms 


we  obtain  finally: 


5 -  rc - - - 


6.  6.  4 
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Note  that  for  frozen  composition,  the  derivatives  of  the  molecular  weight,  IT)  i 
vanish  so  that 

t = i’  <>■ 

f,.  7  Evaliiation  of  Equilibrium  Sound  Velocity 
The  equilibrium  sound  velocity  is  given  by: 

4 


6.7. 


The  above  can  be  expressed  in  terms  of'pjV,  and  Tand  the  specific  heat  data 
by  a  number  of  ways.  One  is  by  the  use  of  the  specific  heat  ratio 

■rl^.  SO  that  eq.(6.  7.  2)  becomes 


Now 


and 


c 


Ye  =* 


6.7,3 


Since,  by  the  partial  derivative  relation  between  three  variables, 

and 

6.7.5 


6.7.4 


Substituting  eq8.(6.  7.  4)and{6.  7.  5)into  eq.  (6.  7.  3)and  simplifying  gives: 


•p  6.  7.  6 


Hence 


Ve^^)r 


Therefore  the  equilibrium  sound  velocity  is  given  by  the  expression 

a* 


6.  7.  7 


6.7.8 
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Substituting  eqs.  (6.  5.  2)  and  (6.  6.  4)  into  (6.  7.  8)  yields  after  rearrangement 


6.7.9 


for  frozen  composition,  =■  [y t}^~  ^  *9*  (6*  9)  reduces  to  the 

classical  thermodynamic  relation 


6.7. 


a"  =  ^r/m 

It  should  be  noted  that  the  equilibrium  velocity  of  sound  is  not  equal  to 
^eRT/  m  .  If  one  wishes  to  use  a  similar  expression  it  can  be  done  by  in¬ 
troducing  a  quantity  P  such  that 


af-VRT/ru 


6.  7.  11 


Then  by  comparing  eqs.  (6.  7.  9),  and  (6.  7.  11) 

1 


r= 


Clearly,  if  the  molar  mass  rn  is  constant 

T”  "  Cp.-V7m  "  W 


6.7.  12 


6.7.  13 


However  if  the  molar  mass  is  varying, 

while  from  the  temperature -entropy  relations,  with  e  the  internal  energy 


Hence 


T= 


which  should  be  contrasted  to  used  in  eq.  (6.  7.  8). 
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7.  COMPUTER  PROGRAM  AND  DATA 

The  program  written  for  the  IBM  7090  digital  computer  at  the  Computer 
Center,  University  of  California,  Berkeley,  to  evaluate  the  detonation  para¬ 
meters  of  gaseous  mixtures  consists  of  four  programs  and  subprograms,  'fhe 
Fortran  listings  are  given  in  Tables  2-7  which  are  discussed  below,  A  flow 
diagram  illustrating  the  operations  performed  by  the  four  sub- routines  and 
their  interdependence  is  shown  in  Fig.  5. 

Table  2  Source  Program 

This  program  loads  the  initial  data  into  the  computer  and  evaluates  the 
upstream  thermodynamic  properties.  These  initial  conditions  are  then  printed 
out,  and  several  constants  are  evaluated  from  the  input  data  which  are  used 
later  in  the  program. 

Tables  3,  4,  5  Subroutine  Spec 

There  are  three  Spec  subroutines  which  shall  be  denoted  by  Spec(a), 
Spec(b),and  Spec(c).  The  purpose  of  the  subroutine  in  all  three  cases  is  either 
to  specify  the  pressure  at  which  properties  are  to  be  evaluated  on  the  Hugoniot 
or  Rankine- Hugoniot  curve,  or  to  converge  to  a  particular  pressure,  e.g.  , 

Pvw  '  O’"  Pv-l- 

The  thermodynamic  properties  and  equilibrium  composition  are  printed 
out  for  each  solution  in  this  subroutine. 

Table  6  Subroutine  Main 

Subroutine  main  is  called  by  the  Spec  subroutine  to  determine  the 
correct  value  of  temperature  on  the  Hugoniot  (or  Rankine -Hugoniot)  curve 
corresponding  to  the  given  value  of  P.  After  converging  to  the  correct  0, 
the  equilibrium  sound  velocity  and  specific  heat  ratio  are  evaluated  and  the 
subroutine  returns  to  subroutine  spec. 
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Table  7  Subroutine  Brink 

This  subroutine  is  called  by  the  subroutine  main  at  each  approximate 
value  of  T and  P  to  evaluate  the  equilibrium  composition  and  several  thermo¬ 
dynamic  properties. 

Table  5  Input  Data 

The  input  data  is  given  to  the  computer  on  IBM  data  cards.  Instructions 
for  putting  the  required  data  on  the  cards  are  given  below.  The  fox  mat  state¬ 
ments  given  by  each  data  card  and  the  numbers  on  the  left  indicate  the 
spaces  on  the  card  where  the  data  must  be  inserted 
Data  1:  1  card  (7E10.  3) 

I- 10:  ^  f  atmospheres 

II- 20:  7/*  degrees  Kelvin 

21-30:  convergence  criterion  for  AZ'J 

31-40:  convergence  criterion  for  ^  ^ j  ;  ^ 

41-50:  increment  of  P  for  Spec  (a) 

blank  for  Spec(b)and  Spec(c) 

51-60:  initial  value  at  P  for  Spec(a) 

initial  approximation  of  for  Spec(b) 

initial  approximation  of  for  Spec(c) 

61-70:  initial  approximation  for  T  corresponding  to  the  initial  P 

Data  2:  1  card  (7£10.  3) 

I- 10:  blank 

II- 20:  initial  approximation  for  for  Spec(b) 

blank  for  Spec(a)  and  Spec(c) 
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Data  3:  1  card  (7110) 

I- 10:  number  of  pressures  to  be  evaluated  on  Hugoniot  curve 

for  Spec  (a) 

blank  for  Spec(b)and  Spec(c) 

II- 20:  number  of  constituents  in  downstream  equilibrium  system 

21-30:  number  of  components  in  downstream  equilibrium  system 

31-40:  number  of  species  in  upstream  mixture 

41-50:  1  if  f^«.|  is  to  be  computed  for  Spec(b) 

0  if  is  not  to  be  computed  for  Spec(b) 

blank  for  Spec(a)and  Spec(c) 

Data  4:  1  or  more  cards  (7110) 

I- 10:  subscripts  of  constituents  appearing  in  Initial  mixture 

Data  5 

II- 20:  subscripts  of  constituents  appearing  in  lists  of  components 

Data  6:  twice  the  number  of  cards  as  the  number  of  constituents  in  downstream 
system 

/ 

1-70:  Data  on  each  card  is  a  line  from  Table  8  :  i.  e. ,  each 

card  has  coefficients  for  one  species  and  for  one 
temperature  range.  The  cards  with  the  coefficients  for 
the  constituents  must  appear  in  the  order  that  the 
constituents  appear  in  the  ,  i.  e. ,  Data  8.  The  low 

temperature  range  coefficients  must  appear  first 
followed  by  the  high  temperature  range  coefficients. 

Data  7:  1  or  more  cards 


10  spaces 
per  data: 


molecular  weights  of  components  in  order 
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Data  8:  1  or  more  cards  (7£10.  3) 

The  matrix  of  the  stoichiometric  coefficients^  ,  with  one  or  more 
cards  for  each  l  .  If  J  becomes  greater  than  7,  use  two  cards  for  each  l  . 
Data  9:  1  or  more  cards 
10  spaces 

per  data:  relative  composition  of  upstream  mixture  in  order  of  Data  4. 
Data  10:  1  or  more  cards  (7A10) 

10  spaces 

per  data:  symbols  of  constituents 
Data  11:  1  or  more  cards  (7£10.  3) 

Approximate  mole  fractions  of  components  in  equilibrium  product^. 

Data  12:  1  card 

blank  for  Spec(a)  and  Spec(b) 

value  of  Mach  number  at  which  shock  wave  properties  are  to  be  computed 
for  Spec(c) 


Table  8  Preparation  of  Thermodynamic  Data 

A  fifth  degree  polynomial  was  fitted  to  the  specific  heat  data  of  each 
species  so  that  the  molar  specific  heats  were  given  by: 

(d  L  ^nCiT-lO  J  /ooi  1^00 

•0*0  ' 

£■  7,6,1 

^  )  isoo  ^sioo 


where  two  series  expression  were  used  to  obtain  a  more  accurate  fit  over  the* 
entire  temperature  range.  The  molar  enthalpy  of  the  species  was  found 
then  by  integration  of  eq,  {7.  6.  1),  yielding  thus: 


Wr"  ^  ciu/t) 


100  Iroo 


7.6.2 

lsooi~r  ^5*160 


i  ‘t'  1  wmlit'  :»•  M 
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where  the  constants  of  integration  and  were  adjusted  so  that  the 
enthalpy  computed  from  both  eqs.  (7.  6.  2)  were  identical  to  the  value  listed 
in  the  JANAF  Ubles  at  1500°K. 

Dividing  eq.  (7.  6.  1)  by  the  temperature  and  integrating  the  resulting 
expression  yielded  finally  an  equation  for  the  molar  entropy: 

3 

where  again  the  constants  and  were  adjusted  to  provide  an  identical 
fit  with  the  JANAF  data  at  1500°K. 

Table  9  lists  three  sets  of  sample  input  data  for  the  computer  program 
using  Spec(a),  Spec(b),  and  Spec(c)  respectively  where  each  set  is  for  a 

mixture  initially  at  1  atm  and  288.  72^K  (60^F).  The  corresponding 
computer  output  is  then  listed  in  Tables  10-  12. 


"1^  =  0.oiloqT  oJt •  ■Co  /ooiTi/Soo 

t/oj7”  T  ^  ^  bwi  (t-  lo'COt  kj- 


7.6. 


S'lOO 
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8.  DISCUSSION  OF  RESULTS 

Detonation  parameters  for  the  forty  different  initial  conditions  listed  on 
page  4  were  calculated  using  the  computer  program  and  the  results  are  pre¬ 
sented  in  Figs.  6-31. 

Figs.  6-9  show  the  P“Vrepresentation  of  the  Hugoniot  curves  computed 
for  initial  pressures  of  .  1  and  760  mm  Hg.  Fig.  6  presents  the  Hugoniot  curves 
for  a  mixture  initially  at  200^F,  Fig.  7  shows  the  Hugoniot  curves 

for  all  three  compositions  at  an  initial*  temperature  of  60°F,  Fig.  8  illustrates 
the  results  for  a  SHz'hOp  mixture  initially  at  -50®F  while,  finally  Fig.  9 
displays  the  Hugoniot  curves  for  all  compositions  at  an  initial  temperature  of 
-180°F.  It  is  of  interest  to  note  that  the  Hugoniot  curves  are  relatively  inde¬ 
pendent  of  initial  composition,  and  further  that  the  CJ  state  in  all  cases  has  a 
volume,  or  velocity,  ratio  of  approximately  0.  54, 

Figs.  10  and  11  show  dimensionless  temperature -entropy  and  pressure- 
entropy  diagrams  for  the  case  of  a  f  mixture  initially  at  1  atm  and 

60®F.  They  illustrate  the  classical  condition  that  the  entropy  is  a  minimum 
at  the  CJ  state.  Further  they  indicate  that  the  pressure  variation  is  greater 
than  that  of  the  temperature  over  a  given  portion  on  the  Hugoniot  curve. 

The  effect  of  initial  pressure  on  the  CJ  pressure  and  temperature  ratios 
are  shown  respectively  in  Figs.  12  and  14.  Fig.  13  is  a  cross  plot  of  Fig.  12 
for  the  composition  and  illustrates  the  obviously  greater  effect  of 

I 

initial  temperature  on  the  CJ  pressure  ratio. 

Figs.  15  and  17  show  the  effect  of  initial  pressure  on  the  CJ  wave  velocity 
and  Mach  number  respectively.  It  is  interesting  to  note  that  the  CJ  wave 
velocity  is  comparatively  independent  of  initial  temperature,  while  the  Mach 
number  is  relatively  independent  of  composition.  Fig.  16»  obtained  from  a 
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cross  plot  of  Fig.  15  for  the  3^2+02  composition  demonstrates  the  effect 

of  initial  temperature  on  wave  velocity. 

Figs.  18-25  show  the  equilibrium  composition  at  the  CJ  state  as  a 

function  of  initial  pressure  for  each  initial  composition  and  temperature. 

They  reveal  the  interesting  result  that  the  equilibrium  composition  of  HjO 

is  approximately  the  same  for  all  conditions. 

Figs.  26  and  27  indicate  the  effect  of  initial  pressure  on  both  frozen 
and  equilibrium  CJ  specific  heats  for  all  three  compositions  at  -ISOOF  and 
60°F.  The  frozen  specific  heat  is  seen  to  be  independent  of  initial  compo¬ 
sition.  More  interesting  yet  is  the  comparison  between  the  equilibrium  and 
frozen  specific  heat,  the  former  being  about  five  times  larger  than  the  latter. 

Fig.  28  illustrates  the  dependence  of  both  frozen  and  equilibrium  CJ 
specific  heat  ratios  on  the  initial  pressure.  The  equilibrium  specific  heat 
ratio,  ,  is  nearly  independent  of  initial  pressure  and  temperature,  but 
shows  a  significant  dependence  on  composition.  The  frozen  composition 
specific  heat  ratio.  .  on  the  other  hand,  shows  an  obvious  dependence 
bn  initial  pressure,  although  it  is  relatively  independent  of  initial  tempera- 

ture  and  composition. 

The  von  Neumann  spike  values  of  P  ,  ©  ,  and  V  are  shown  in  Figs.  29. 
30,  and  31  as  a  function  of  initial  pressure.  It  is  interesting  to  note  that  the 
ratio  of  the  VN  pressure  and  temperature  to  the  corresponding  CJ  values  is 
nearly  independent  of  the  initial  conditions. 

The  CJ  properties  for  the  initial  conditions  considered  are  tabulated 
in  Table  13,  while  the  von  Neumann  spike  properties  are  presented  in 

Table  14. 

Finally,  it  is  of  interest  to  compare  results  of  the  present  calculations 
vvith  those  obtained  by  other  investigators  together  with  some  experimental 
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results.  The  data  for  CJ  detonations  presented  in  Table  15  indicate  a  good 
agreement  between  the  various  results  based  on  the  use  of  the  equilibrium 
sound  speed,  which  in  general  are  slightly  larger  than  both  experimental 
measurements  and  those  results  based  on  frozen  sound  speed.  The  present 
results  also  agree  quite  well  with  those  of  Bollinger  and  Edse  (7),  which  were 
obtained  at  313®K. 
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TABLE  1 


Comparison  of  Various  Methods  of  Analysis 


INDEPENDENT 

author  system  scope  variables  method  remarks 


1 - 

Lewis  & 
Frleuf  (1) 

Hydrogen-Oxygen 
with  Inerts 

CJ  state 

lE^li 

Iteration  onT2  ,Vp  ,  to 
satisfy  conservatism  equa¬ 

tions  with  CJ  state  deter- 
mlni^  fi^  cradltlon 

- 

Atomic  oxygen  not 
considered  In  the 
eqxilllbrlum  system 

Bsrsts , 
et  «1  (2) 

Hydrogen , Oxygen 

CJ  state 

Same  as  Lewis  and 

Frlauf 

Edse  (3) 

Arbitrary 

Hugonlot 
curve  and 
CJ  state 

^2 

Hugonlot  curve  determined 

1^  selecting  arbitrary 

J2  end  Iterating  for  the 
corresponding  value  o£p2  • 

CJ  state  determined  whih 
slope  of  Hugonlot  curve  Is 
Identical  to  slope  of 
Rayleigh  line. 

Wolf son 
&  DuTtn  (4) 

Arbitrary 

CJ  state 

Tg.Pz 

Iterations  on  \  P2 

satisfy  conservation  equa¬ 
tion  with  CJ  point  determin¬ 
ed  from  condltlonM2-M2^1<0 

Ltjksr 
•t  al  (5) 

Hydrogen, Oxygen 
and  Steam 

CJ  state 

T2.P2 

Same  method  as  Reference  1 

Van  de  Waal  equa¬ 
tions  of  state 
used  for  H2O 

Elscn 
at  al  (6) 

C,H,0,N,A 

■ 

Ml 

For  selected  values  of  , 
Iterations  are  preformed 
to  determine  points  along 
Hugonlot  curve.  CJ  state 
specified  by  condition  that 
M<|  Is  a  minimum. 

First  to  check  CJ 
state  %nth  Mp^'l 

Bollinger 
&  Edse  (7) 

Hydrogen-Oxygen 

CJ  state 

^2 

Hugonlot  curve  foimd  by 
selecting  arbitrary  To  ,  and 
determining  corresponding  Pp 
by  Iteration.  The  CJ  state 
Is  determined  by  finding  the 
value  of  minimum  . 

Equilibrium  equa¬ 
tions  are  Incorpor¬ 
ated  directly  In  the 
conservation  equa¬ 
tions 

Barrera 
et  al  (8) 

Arbitrary 

CJ  state 

P2'^2 

Computes  equilibrium  compo¬ 
sition  for  several  values 
of  p  and  T  In  the  neighbor¬ 
hood  of  the  CJ  state.  Values 
of  Tq  and  h  are  then  evalu¬ 
ated  from  both  the  conserva¬ 
tion  equations  and  the  equa¬ 
tion  of  state.  A  cross  plot 
technique  yields  the  solu¬ 
tion  to  the  CJ  state.  The 
criteria  for  the  CJ  state  Is 
that  of  minimum  entropy. 

Graphical  method 

Zeleznlk  & 
Gordon  (9) 

Arbitrary 

CJ  state 

P2.T2 

Iteration  on  pu  and  Ip  vlth 

condition  M2=  ^Mpe  ^corpor- 

ated  In  confervaribn 
eqxiatlons 

Sound  speed  of 
proci^ts  expressed 
by  c^=  TRT  ^  where 

'  IsKw  Cje 

Bird 

et  al  (10) 

Arbitrary  Includ¬ 
ing  condensed 
phase 

Hugonlot 
curve  & 

CJ  state 

^2 

^  Is  determined  Itera¬ 
tively  for  specified  T2  to 
satisfy  Hugonlot  eqxiatlon 
and  equation  of  state.  CJ 
date  evaluated  for  M<| 
minumum. 

May  Inclixle  more 
than  one  phase  In 
the  chemical  system 

Present 

Work 

Arbitrary 

CJ  state, 
von  Neumann 
spike  and 
Hugonlot 
curve 

P2 

Tp  Is  determined  Itera¬ 
tively  for  specified  P2  to 
satisfy  Hugonlot  equation 
and  equation  of  state.  CJ 
state  evaluated  forM^2’^* 

Sound  speed  of 

V  .  Cae 

TABLE  2,. 


Fortran  Listing  -  Source  Program 


*  LIST 

*  LABEL 

*  roRfirWt^ -  ^ 

DIMENSION  Aa0tl0)»B(20tl0)tBe(20)»BC(20*10)»CI30«10)«COMP(20) t 

1CPM(20)»AEQ(10*10)>OZ(20)*G(20)*G6(10) *H(20) tHF(20) *0( 10) tUU(20) . 
2UB(10)»WMM(20)*NB(10)*NX(20)*XI20>*Z(10)*WMC(10)*WMP(20)*CP(20)t 
3ST0RE ( 50 120 ) (BETA ( 20 ) tDH( 20 ) 

COMMON  AtAA*AAAtB*BB«BBB«BC*C«CC*C0MP*C0Z*C0V*CPM.CPl*CP2*0*DP»0Tt 

-  10T2*DZ*6*GAMMAl»6&tH*HF»Hl»H2iMM*MP»NB«NCM»NC0M>NC0NtN60tNN*NX*P> 

2Pl>P2.0t00*R*S0NIC>T.Tl»T2*Ul*UU«UB*VtVltV2.WMtWMM»WM0»X»XN.Z.ZZ. 

3E*WMl»WMC*WMP*CONV*CON*CP*EMACH2*FMACH2*XMCHl*FRSON«EOSON*GAMM2Rt 

4GAMM2EtCPR*CPE»H2*ST0RE*U2*CVE*CVR*BETA»0H*LSMFT*DStDU.0ELP«XX»YY> 

5LX.LY»LZ*LW*LV*LU»LT»LS*LR*LD*0L0SStS0UND»CHECKV»CHi:CK*N»RR*RRM.SS 

1000  FORMAT  (7E10.6) 

1001  FORMAT  (7110) 

1003  FORMAT  (7A10) 

1004  FORMAT  (6E12.5) 

1005  F0RMAT(8E10.3) 

R>1. 98726 
C0NV=4184. 

-  . -  C0N*41#84/L*013» 

1050  READ  10O0*PltTl*CDZ«COV.OP*P*TtQLOSS*CHECK 
NGO«0 
NN^O 

IF(Pl)  1100*2000*1055 
1100  PRINT  1110 

- mo  FORMAT  «7H  STOP  1)-  - 

CALL  EXIT 
1055  PRINT  1060 
ARG6S0 

1060  F0RMAT(56H1HUG0NI0T  CURVE  CALCULATIONS  FOR  GASEOUS  MIXTURES - PROP 

162HULSION  dynamics  LABORATORY— UNIVERSITY  OF  CALIFORNIA— BERKELEY) 
— 1150-ReAD-100t»MP^rN€OH»NCOH,NCR*LR 
READ  1001*  (NX(I)*l3l*NCM) 

READ  1001*  (NB(J)*J>l*NCOM) 

N>NCON 

DO  1160  Isl*NCM 
IF  (N-NX(I))1155*1160*1160 

- — tt58  7**NXM)-  - 

1160  CONTINUE 
MC»2«N 

READ  1005*((C(I*J)*J»1*8)*I*1*MC) 

1200  READ  1000*  (WMC( J) * J«l*NCOM) 

DO  1250  I>1*N 

- rtSfr-READ  10O0*tB(+*;l)*J»l*NC0M) 

READ  1000*(X(I)*I«1*NCM) 

READ  1003*(COMP(J)*J>1*N) 

00>0*0 

DO  1300  I«1*NCM 
1300  0Q-OO4X(I) 

-  DO  1350  I*1*NCM 
1350  X(I)>X(I)/QO 

OQ>0.0 

DO  1500  J>l*NCOM 
0( J)*0*0 


TABLE  2  ~  Continued 


IsNX(K> 

1450  Q( J)»Q< 1 *J>*X(K) 

1500  QQsQQ^QjJ)  _  _  _  _  ,  ... 

00 '1550  JaltNCOM 

1550  0«J)*Q«J1/QQ 
WM0«0 

00  1560  JsltNCOM 
1560  WMO=WMO-*-WMC«  J)«0<  J) 

WM1*WM0»00 
■  "RRsR/WMi 
RRM=RR*T1 
'00  1690  Isl.NCON 
WMP( 1 )»0* 

00  1690  J«1.NC0M 
1690  WMP( 1 )»WMP< 1 )+B( 1 tJ»*WMC( J» 

00  1700'  1  =  1.NC0N 
6B( I )»0.0 
DO  1700  J=ltNCOM 
1700  BBC  I  )»BB«  n+B( !  t  J> 

DO  1750  IsltNCON 
DO  1750  J=1*NC0M 

'  1750'  BC( I  * J)=BC I *J)-CBBC I )-l. )*Q< JI 
TK»T1/1000. 

1760  IFCTl-1500.)  1770*1770*1772 

1770  00  1771  K=1*NCM 

1771  CPmMM  C  { ( CC(  I  *6)»TX+C(  1  *5»  t*TX-*-CC  I  *4)  )*TK+C(  I*3))*TX'C-CCI*2>»*TK+ 

1C(1*1))*R' 

GO  TO  1780 

1772  DO  1773  K=1*NCM 
IsNXCIO 

1773  CPmI!c)  =  (  ( C  ( (  C(J*6)*TX'fCC  J*5)  l^TK+CC  J*4)  (•TX-t-C C  J* 3 )  )*TX+C  ( J* 2 ) ) »TIC+ 

lT(J*n»»R 


1780  CP1«0. 

DO  1800  I»1*NCM 
1800  CP1«CP1+CPM( I )*X( 1 » 

GAMMA1*CP1/CCP1-R) 

CP1»CP1/WM1 

"  ■■■  S0NIC=GAMMA1'*R*T1/WM1 
U1«S0RTFCS0N1C*C0NV) 

V1*R*T1/(WM1*P1 ) 

V*V1*C0N 

1F(T1-1500»)  1810*1810*1812 

1810  DO  1811  ^ 

HFc'xj!^(C(((C<l*6J*TX/5.+CCI*5»/4.l»TKKCl*4)/3.)*T<fCCI*3)/2.)*TX+ 

ICC  I *2  > l^TK+C  C I  *  1 )4LOGF  C  T1 ) C 1 *8 )*R 

1811  HCX)  =  CCC<CCCI  *6  J^LW^/B^+CC  I  *5)/5*  )*TX+CC  I  *4)  /4*  )*TX*CC  I  *3 ) /3* 

ICC  I  *2  J/2.  )*TX*'CC  1*1)  )*T1*R+CC  1  *7)*R 

60  TO  1860 

"1'812  OO  1813  K«1*NCM 
I«NXCX) 

HFCK)sC  C  C  C  CCCJ*6)*TX/5*+CCJ»5)/4*  j^TX+C  C  J*4 )/3* ) *TX*C  C J*3 ) /2*  *TX  + 
1CCJ*2))*TXKCJ*1)«L06FCT1)  )*R+CCJ*8)*R 

1813  HCX)»CCCCCCCJ*6)*TX/6.+CCJ*5)/5.)*TX+CCJ*4)/4.)*TX4CCJ*3)/3.I*TX4 

- jV2)72«)*TX*CI  J*l)  )*TI»R+CC  J»7)»R 

1850  H1«0. 

SS>0* 


TA.  Lli  2  -  Continued 


00  1900  J«1»NCM 

SS«SS+  X(J»*<HF<J»-R*L0GF(P1*X(J)) ) 

H1<H1/WM1 
SS«SS/R 
HH«H1/RRM 
HC=CP1/RR 
PRINT  1950 

iW-F0RMATT50H0  HHTiAL  CONDITION  OF  MIXTURE) 

PRINT  1960fPl*Tl«HH.Ul 

1960  format (4H0PlsE12«5*  14HATM  TEMP«F10.4*16H  DECK 

/1.2«16H  SO  VEL=F10.4*5HM/SEC) 

PRINT  1961 (V  *  HCtGAMMAltWMl 

1961  F0RMAT(4H  V1>E10*5*16HCM3/GM  CP1-F10.3*16H 

l#3fl6H  WM1-F10.3/) 

“RINT  1959»SS 
1959  FORMAT (9H  ENTROPY»F10.4 ) 

PRINT  1962 

1962  FORMAT(23HOCOMPOS1TION  OF  MIXTURE/) 

DO  1965  J-l.NCM 

- . i-NX+U^)- - 

PRINT  1963«COMP(ll*X(J) 

1963  F0RMAT(A7*2H  »E12«5) 

1965  CONTINUE 

READ  1000*  (X(J)*J«l*NCOM) 

1980  CONTINUE 

- -CAtt-SPEC- -  ...  - 

GO  TO  1050 
2000  CAL(  EXIT 
END 


ENTHALPY*F10 

GAMMAlsFlO 


TABLE  3.  Fortran  Listing  -  Spec(a)  Subix>utlne; 
Hugonlot  Curve  Calctilations 


*  LIST 

*  LABEL 

*  FORTRAN 

_ SUBROUTINE  SPEC _ _ _ _ _ 

DIMENSION  A( 10»10 I •R( 20» 10 ) «BB( 20 ) «BC( 20»10 > tC ( 30«10) tCOMPC 20 ) • 
1CPM(20).AEO(10.10).OZ(20).G(20).GG(10).HI20).HF(20).QI10)*UU(20)» 
2UB(10)«WMM(20)»NB(10>«NX(20)>X(20)«Z(10)«WMC(10>»WMP(20)>CP(20)» 
3ST0RE(50.20) .BETA(20) ♦0HI20) 

COMMON  A»AAtAAA»B»BB»BBB»BC»C»CC*C0MP.C0Z,C0V.CPM,CPl.CP2»D.DP.0T» 
1DT2.DZ.G»GAMMA1»GG,H»HF»H1,H2»MM»MP»NB%NCM,NC0M»NC0N»NG0»NN»NX»P» 
2Pl.P2.0.QQ.R»S0NIC.T.Tl.T2»Ul,UU.UB,V.Vl»V2.WM,WMM,WM0.XtXN»ZtZZ» 
3E . WMl ,WMC .WMP.CONV .CON .CP .EMACH2 . FMACH2 . XMCHl . FRSON • EQSON .6AMM2R • 
4GAMM2E»CPR.CPE»H2»ST0RE»U2.CVE.CVR»BETA.DH.LSMFT.0S.0U.DELP,XX.Yy. 
5LX.LY, LZ.LW.LV.LU.LT.LS.LR.LD.OLOSS. SOUND. CHECKV.CHECX.N.RR.RRM.SS 

1  FORMAT (IH  ) 

2  F0RMaT(47H1THFRM0DYNAMIC  PROPERTIES  ON  THE  HUGONIOT  CURVE) _ 

3  format (56H0  PRESS  R  TEMP  R  VEL  R  MACHl  FMACH2  FMACH2  FSDVEL2 

165HESDVEL2  FGAMM2  EGAMM2  FSPHT  ESPHT  ENIH2  MOLWT  ENTROP 
2Y  ) 

4  FORMAT(F9.4.F8.4.F7,4.3F7,4.2F9.3.2F7.4.2F8.4.F10.3.F8.3.F9.4) 

5  FORMAT(F10.4.9E12.5) 

ft  F0RMATt46HlEQUILIBRtUM  COMPOSITION  ON  THE  HUGONIOT  CURVE) _ 

7  FORMAT! lOHO  PRESS  R  9A12) 

NGO=0 
2000  MM=5 

LSMFT=0 

K=NCON+l 

LV=0 

- PftiNt  5 - 

PRINT  1 
PRINT  3 
GO  TO  2250 
2200  P=P+0P 
2250  P2=P*P1 

- 12*1*11 - 

CALL  main 

XX=XX/RR 

H2=H2/RRM 

CPE=CPE/RR 

CPR=CPR/RR 

- IFIHW-5)  2300.2275 '72T5 - - 

2275  PRINT  1 
MM=0 

2300  PRINT  4. P.T.V. XMCHl. FMACH2.EMACH2. FRSON. EOSON.GAMM2R.GAMM2E. 
1CPR.CPE.H2.WM,XX 
MM=MM+1 

- NGOsNGO+l - — - - 

LSMFT=LSMFT+i 
STORE (LSMFT.1)=P 
DO  2700  I=l.NCON 
L=I  +  1 

2700  STORE(LSMFT.L)=X(  I  ) 

- IF(MP-NOe)  i000>2800l2750 - - 

2750  IF(40-LSMFT)  2800.2800.2200 
2800  MM=5 

PRINT  6 


TABLE  a 


Fortran  Listing  Spec(b)  Subroutine;  CJ  Calculations 


*  LIST 

*  LABEL 

*  FORTRAN 
SUBROUTINE  SPEC 

DIMENSION  A ( 10,10) ,6(20*10) tBB( 20) *BC( 20*10) *C( 30*10) *C0MP(20) * 
1CPM(20)*AEO( 10*10) *DZ(20) *G(20) *GG( 10) *H(20) *HF(20) *0? 10) *UU(20) * 
2UB( 10)*WMM(20) *N3( 10) *NX(20) *X(20) *Z ( 10 ) *WMC ( 10) *WMP( 20 ) *CP( 20 ) * 
3STORE(S0*20) *BETA ( 20 ) *0H ( 20 ) 

COMMON  A*AA*AAA*B*BB*BBB*BC*C*CC*C0MP*CDZ*CDV*CPM*CP1*CP2*D*DP*DT* 
10T2*DZ*G*GAMMA1 *6G*H*HF*H1 ,H2 *MM *MP *Nb *NCM ,NC0M * NCON *NGO*NN *NX *P , 
2P1*P2*0*QQ*R*S0N1C*T*T1*T2*U1*UU*UB*V*V1*V2*WM,WMM*WM0*X*XN*Z*ZZ* 
3E*WM1,WMC*WMP*CONV*CON*CP*EMACH2*FMACH2*XMCH1*FRSON*EQSON*GAMM2R* 
4GAMM2E*CPR*CPE*H2*ST0RE*U2*CyE*CVR*BETA*DH*LSMFT*DS*DU*DELP*XX*YY* 
5LX*LY*LZ*LW*LV*LU*LT*LS*LR,LD*OLOSS*SOUNO*CHECXV*CHECK,N*RR*RRM*SS 
IF  (LR)  5*5*4 

4  PX=0. 

GO  TO  2100 

5  PP  =  P 
TT  =  T 
P=CHECK 
P2=P*P1 
T2=T*T1 
LV--=1 

10  CALL  MAIN 
DS=H2 
DD=H2 

BBBB=HUV1#(P2-P1 1-rULOSS 

PsP+,05 

P2=P#P1 

CALL  MAIN 

DS*(H2-DS) 

AAA=V1»P1 

CC*DS/*05 

0ELP=(DD-BBBB)/(AAA-CC) 

IF(ABSF(DELP)-1.»  14*12*12 
12  0ELP=0ELP/ABSF(DELP) 

14  CONTINUE 

P=P-,05+DELP 

I F ( ABSF ( BBBB-OD ) -CDV ) 20  *  20  *  1 0 

20  IF( ABSF(DELP)-CDV)  21*21*10 

21  PX  =  P 

PRINT  22*  PX 

22  F0RMAT(63H0  THE  PRESSURE  RATIO  OF  THE  HUGONIOT  CURVE  AT  V2=V1  IS  E 
IQUAL  T0F9r4) 

IF(PX-1.)  23*23*28 

23  PRINT  25 

25  FORMAT (48H0THERE  IS  NO  CHAPMAN- JOUGUET  POINT  FOR  THIS  DATA) 

GO  TO  3270 
28  P=PP 
T  =  TT 

DO  2000  1  =  1, NCM 
K=NX( I ) 

2000  DH( I )=X(K) 

00  2050  1=1, NCM 
2050  X( I )=DH( I ) 

2100  PRINT  2200 
2200  FORMATdH  //) 

PRINT  2250 

2250  F0RMAT(54H  THERMODYNAMIC  PROPERTIES  AT  THE  CHAPMAN- JOUGUET  POINT) 
P2=P*P1 
T2=T*T1 


TABLE  4  -  Continued 


LV«0 

2300  CALL  MAIN 

- ^«EMACH2 -  -  -  - 

PP*P 
P»P**-«05 
P2»P»P1 
SOUND*EMACH2 
CALL  MAIN 

—  _ _ _<VH  «■  t  /V|»  ._.... 

t7U*Tc::fiACrnr*t7t7xA«t7!> 

t>ELP=(l. -SOUND! /DU 
IF<ABSFCDELP)-1.)  2310t230e*2308 
2308  DELP=DELP/ABSFJDELP) 

2310  CO((lTINUE 
P=PP+DELP 
-  P2*P*P1 

IF<ABSF< SOUND-1. )  -CDVI  2350*2350*2300 
2350  IF(ABSF(OELP)-CDV)  2500*2500*2300 
2500  CONTINUE 
CALL  MAIN 
T=T2/T1 
PRINT  2520 

2520  FORMAT (56H0  PRESS  R  TEMP  R  VEL  R  MACHl  FMACH2  EMACH2  FSDVEL2 
165HESDVEL2  FGAMM2  EGAMM2  FSPHT  ESPHT  ENTH2  MOLWT  ENTROP 
2Y  ) 

PRINT  2560 
XX=XX/RR 

- - H2-H2/RRM 

CPE»CPE/RR 

CPR*CPR/RR 

PRINT  2550*P*T*V*XMCH1*FMACH2*EMACH2*FRS0N*EQS0N»GAMM2R»6AMM2E* 
1CPR*CPE*H2*WM*XX 

2550  F0RMAT(F9.4*F8.4*F7.4*3F7,4*2F9.3*2F7.4*2F8.4*F10.3*F8.3*F9.4) 

- 25*0- FORMAT IlH  I 

PRINT  2560 
PRINT  2600 

2600  F0RMAT<54H0  EQUILIBRIUM  COMPOSITION  AT  THE  CHAPMAN- JOUGUET  POINT) 
PRINT  2560 
LFIRST=1 

IF<NCON-9)  3000*3000*2950 
2950  LPRINT=9 

GO  TO  3050 
3000  LPRINT=NCON 

3050  PRINT  3060*(COMP( J)*J=  LF IRST *LPRINT ) 

3060  F0RMAT<9A12) 

-  .  PRINT  2560 

PRINT  3070* <X( I ) *IsLFIRST.LPRINT) 

IF(NCON-LPRINT)  3270*3270*3260 
3070  F0RMAT(9E12.5) 

3260  LFIRST=LFIRST+9 

IF(NC0N-LPRINT-9)  3262*3262*3264 

^  _  - ^42  LPRINT»NC0N 

GO  TO  3266 

3264  LPRINT=LPRINT+9 
3266  PRINT  2560 
GO  TO  3050 
3270  RETURN 

'■ . .  END 


TABLE  3»  Fortran  Listing  -  Spec  (c)  Subroutine,  Shock  Conditions 
at  Given  Mach  Number 


*  LIST 

*  LABEL 

*  FORTRAN 

- SUBRQUI1!IE_SPK _ _ 

DIMENSION  A( lOtlOl *B(20«10 I *BB(20) *BCI 20*10) *C( 30*10) *COHP( 20)* 
lCPM(20i*AEO(10*10)*02l20)*6(20)*6G(10) *HI20) *HF( 20) *Q( 10) *UU(20) * 
2UBaO)*WHM(20)*NB(10)*NX(20)*XI20)*Z(10)*WMCI10)*WMPl20)*CPI20)* 
3ST0REI30*20)*BETA(20)*DH(20) 

COMMON  A*AA*AAA*B*BB*BBB*BC*C*CC*C0MP*C0Z*CDV*CPM*CP1*CP2*D*DP*0T* 
10T2*DZ*6*GAMMA1*GG*H*HF*H1*H2*MM*MP*NB*NCM*NC0M*NC0N*NG0*NN*NX*P* 
2P1*P2*Q*00*R*s6nIC*T*T1*T2*U1*UU*UB*V*V1*V2*WM*WMM*WM0*X*XN*Z*ZZ* 
3E*WM1*WMC*WMP*CONV*CON*CP*EMACH2*FMACH2*XMCH1*FRSON*EOSON*GAMM2R* 
4GAMM2E*CPR*CPE*H2*ST0RE*U2*CVE*CVR*BETA*DH*LSMFT*DS*DU*DELP*XX*YY* 
eLX*LY*LZ*LW*LV*LU*LT*LS*LR*LO*OLOSS*SOUNO*CHECKV*CHECK*N*RR*RRM*SS 
1  F0RMAT(8E10*3) 

J^;^0 

. . READ  i*XMCHl 

ALPHA«GAMMA 1«XMCH 1**2 
T2«T*T1 
9000  P2«P*P1 

CALL  MAIN 

__  _ VRHM/ _ 

vTlP)^*  1 .-( P-l .) /ALPHA 

P«P+*05 

P2«P*P1 

CALL  MAIN 

DVDPsIV-VRH)/.05 

_ _ _ DELP« ( VALPH-VRH ) / ( DVDP»1 . /ALPHA ) 

IF(ABSF(DELP)-i«)  9106*9100*9050 
9050  OELP>OELP/ABSF(DELP) 

9100  P»P-.05+0ELP 

I F ( ABSF  <  DELP ) -COV )  9200*9200*9000 
9200  P2*P»P1 

_  CALL  MAIN  _ 

PRINT  9400 

9400  FORMAT (29HOPROPERT1ES  BEHIND  SHOCK  WAVE) 

PRINT  9410 

9410  FORMAT (56HO  PRESS  R  TEMP  R  VEL  R  MACHl  FMACH2  EMACH2  FSDVELl 
165HESDVEL2  FGAMM2  EGAMM2  FSPHT  ESPHT  ENTH2  MOLWT  ENTROP 

_ ZY.J _ _ _ _ 

PRINT  9430 

PRINT  9420*P*T*V*XMCH1*FMACH2*EMACH2*FRS0N*E0S0N*GAMM2R*GAMM2E* 
1CPR*CPE*H2*WM*XX 

9420  FORMAT (F9*4*F8*4*F7*4*3F7*4*2F9*3*2F7*4*2F8*4*F10*3*F8.3*F9*4) 
PRINT  9430 

_ 9>30. FORMAT (IH  ) 

PRINT  9440 

9440  FORMAT (30HOCOMPOSIT ION  BEHIND  SHOCK  WAVE//) 

LFIRST*! 

IF(NCON-9*  9500*9500*9450 
9450  LPRINT*9 
60  TO  9550 
*9500  LPRINT-NCON 

9550  PRINT  9560*(COMP( J) *J=LF1RST*LPR1NT> 

9560  FORMAT(9A12) 

PRINT  9430 

PRINT  9570*(X(I ).IsLFIRST*LPRINT) 

9570  F0RMAT(9E12.5) 

.  fFINCON-LPRINT)  9770*9770*9760 

9760  LFIRST*LFIRST+9 


TABLiE  5  -  Continued 


IF(NC0N-LPRINT-9)  9762.9762 f9764 
9762  LPRINTsNCON 
GO  TO  9766 

9764  LPRINT*LPRINT49 
9766  PRINT  9768 
9768  FORMAT! IN  //) 

GO  TO  4550 
9770  RFTURN 
END 


TABLE  6.  Fortran  Listing 


Main  Subroutine 


*  LIST 

*  LABEL 

*  FORTRAN 
SUBROUTINE  MAIN 

DIMENSION  A(  10*10)  ♦Bi  20*10 )  fB0<  20)  *bC  (  20 *  10  )  *C  ('30 *  10 )  *C0MP  (  20 )  * 
1CPM(20) *AEQ( 10*10) *DZ(20) *0(20)*0G( 10) *H(20) *HF(20) *Q( 10) *UU(20) * 
2UB( 10) *WMM(20) *NB (10 ) *NX ( 20 ) *X ( 20 ) *Z ( 10) *WMC( lO) *WKP(20) *CP(20) * 
3STORE ( 50*20) *BfcrA( 20) *0H( 20) 

COMMON  A*AA*AAA*D*cJB.Bbb*BC*C*CC*COMP*CL»2*CDV*CPM*CPl*CP2*D*DP*DT* 
1DT2*DZ*G*GAMMA1 *GG*H*HF*H1 *H2*MM*MP*Nb*NCM*NCUM*NCuN*r4GO*NN*NX*F* 
2P1 *P2*0*00*R*50NIC*T*  ri *12*01 *UU*UB*V* VI *V2*WM*WMM*WM0*X*XN*Z*ZZ* 
3E*WM1 *WMC*WMP*C0NV*C0N*CP*EMACH2*FMACH2*XMCH1*FRS0N*EQS0N*GAMM2R* 
4GAMM2E  *CPR *CPE *H2 *  STORE *U2 *CVE *CVR *BE T A *DH *L5MFT *0S *DU *DELP *XX * YY * 
5LX*LY*LZ*LW*LV*LU*L1  *LS*LR*LO*OLOSS*SOUNi>*Ch2C>;V*CHcCK*:m*RR*RRM*SS 
4000  CALL  BRINK 
PVTP=0. 

DO  4005  I=1*NC0N 

4005  PVTP  =  PVTP+WMP( I )*X( I  )*bFTA(  I  ) 

PVTP= ( 1,-PVTP/WM)*R/ (WM*P2  ) 

AAA=2.*CP2/(V1*(P2-P1) ) 

,  BBH=2.^(H2-Hl-OLOSS )/( ( P2-P1 ) *V 1 ) -1 • 

CC=PVTP/V1 
D  =  R*T2/(V1*P2<^WM) 

DT2=(BBB-i))/(CC-AAA) 

IF( ABSF(DT2)-400, )  a008 *4008 *400 7 

4007  0T2=400**(DT2/(ABSF(DT2 ) ) ) 

4008  T2=T2+DT? 

IF  ( ABSF(BBB-D)-CDV)  4^10*4000*4000 
4010  DT2=DT2/T1 

IF( AR5F(DT2 )-CDV )  4150*4000*4000 
4150  V=BbB 

T=T2/T1 

IF(LV)  4200*4200*4500  ‘ 

4200  U1  =  2«*(H2-H1-QL0SS)  /  ( l•-V^^*2) 

U2=(VM2)*U1 

XMCHl =SQRTF ( Ul/SONIC ) 

GAMM2R  =  CPR/(CPR-R/W'^) 

F  R  S  0  N = G  A  M :  ^  2  R  *  K  *  T  2  /  W  M 
FMACH2=sbRTF(U2/FRSON) 

FRSON^SQRTF ( FRSON*CONV ) 

DO  4300  J=1*NC0M 
4300  GG(J)=-0(J) 

ARG6=0 

LLL=ISIMEQ( 10*NCOM*1*A*GG*ARG6*OZ) 

IF(LLL-l)  4320*4320*4310 

4310  PRINT  4311 

4311  F0RMAT(7H  STOP  5) 

CALL  EXIT 

4320  DO  4350  I=l*NCON 
BETA( I)=0. 

DO  4340  J=l*NCOM 

4340  BETA( I )=BETA( I ) -B ( I  * J ) *4 ( J , 1 ) 

4350  BETA( I )=BETA( I )-l. 

PVPT=0* 

DO  4400  I=l*NCON 

4400  PVPT  =  PVPT^*WMP(  I  )#X(  I  )»BETA(  I  ) 


TABLE  6  -  Continued 


PVPT  =  ( i.+PVPT/WM) 

gamm?e=pvtp-*2»wm*p?*<^;>/  ( R*p\/Pi*cPt ) 
GAM.y2E=l./(  l«-GAMK2t  ) 

E Q SON  =  G AMM 2 F. ♦R T 2 /  (  P VH T  W '<) 

EMACh2  =  5QRTF  (U2/EvJbUN) 

EQSON=SORTr ( EOSON*CONV ) 

XX  =  0. 

DO  4430  I=l*NCON 

IFJX(  I  )*10.**10-1*  )  44‘j0»44DU*4430 
44  30  XX  =  XX4X(  I  )»(Hr  (  I  )-K^rL0GF  (X(  I  )  ) 

4430  CONTINUE 
XX=XX/WM 
4300  RETURN 
END 


TABLE  7 


Fortran  Listing  -  Brink  Subroutine 


*  LIST 

*  LABEL 

*  FORTRAN 
SUBROUTINE  BRINK 

DIMENSION  A( lOtlO) •B(20«10)«BB(20) •BC(20tl0) •C(30tlO)*COMP(20) • 
1CPM(20) fAEOI 10.10) .02(20) .GI20)»CG(10)»H(20)fHF(20) tO(lO) •UU(20»t 
2UB(10).WMM(20).N9(10).NX(2O).X(20).2tl0).WMC(l0).WMP(20).CP(20). 
3ST0RE(50.2C) .BETA(20).DH(20) 

COMMON  A.AA.AAA.B.BB.BBB.BC.C.CC.COMP.CD2.CDV.CPM.CP1.CP2.D.DP.OT. 
10T2.D2.0.GAMMAl.G0.H.MF.Hl.H2.MM.MP.Nb.i'.CM.NC0M.NC0N.N00.NN.NX.P. 
2Pl.P2.U.(J(j.R.SONlC.T.Tl.T2.Ul.UU.UB.V.Vl.V2.MM.WMM.MMO.X.XN.2.ZZ. 
3E.WM1.WMC.WMP.CONV.CON.CP.EMACH2.FMACH2.XMCH1.FRSON.EOSON.GAMM2R. 
4GAMM2E. CPR. CPE. H2. STORE. U2. eve. CVR. BETA. DH.LSMFT.DS. DU. DEEP. XX. YY. 
5LX. LY.LZ.LW.LV.LU.LT.LS.LR.LD.QLOSS. SOUND. CHECKV. CHECK. N.RR.RRM.SS 
L-0 

TK*T2/1000. 

IF(T2-1500.)  500.500.502 

500  DO  501  I=l.NCON 

CPI n=( ( ( < (C( I.6)*TK+C( 1.5) I*TK+C( I .4) )*TK+C( 1 .3 ) )*TK+C ( I .2 ) )*TK+ 
1C(I.1))*R 

H(I)s(((l(C(I .6)*TK/6.+C( I .5)/5. )*TK+C( 1 .4 ) /4. ) ♦TK+C ( I . 3 ) /3» )*TK+ 
1C(1.2)/2.)*TK+C( I .1) )*T2*R+C( I.7)*R 
HF(I)=(((((C(I.6)*TK/5.+C(I.5)/4.)*TK+C(I.4)/3.)»TK+C(I.3)/2.)*TK+ 
1C( 1 .2  )  )*TK+C( I.l )*L0GF( T2) )*k+C( I .8)*R 

501  UU( I )=H( I )-T2*HF( I ) 

GO  TO  504 

502  DO  503  I=1.NC0N 
J*I+N 

CPI  I  )  =  (  ( (  (  (C(J.6)*TK+CI  J.5)  )»TK+C(  J.4)  )  *7K+C  (J.  3  )  ji«TK+C  IJ  .2 )  )*TK+ 
1C(J.1))*R 

HI  I  )  =  I  I  1 1 ICI J.6)*TK/6.+C(J.5)/5.)*TK+C( J.4)/4.)*TK+C(J.3)/3.)*TK+ 
ICI J.2)/2.)*TK+C(J.l) )*R*T2+C(  J.7)*R 
HFI I  )  =  ( I  I  I (CIJ.6)*TK/5.+Cl J.5)/4.)*TK+C(J.4)/3.)*TK+C(J.3)/2.)*TK+ 
ICS  J.2)  )*TK+CIJ.1)*L0GFIT2)  )*R-i-CI  J.8)*R 

503  UUI I )=H( I )-T2*HF( I  ) 

504  DO  505  J=1.NC0M 
I=NBI J) 

505  UBIJ)=UU(I) 

DO  520  I=1.NC0N 
GI 1 )=0.0 
DO  510  J=1.NC0M 

510  GI 1 )=B( I .J)*UB( J)/IR*T2 )+G( 1 ) 

520  GI I )=G( I )-LOGF (P2)-UUI I )/ (R*12) 

IFINN)  523.523.521 

521  DO  522  J=1.NC0M 

522  X(J)=WMM(J) 

523  DO  531  J=1.NC0M 

531  Z( J)=-L0GFIX(J) »-L0GF(P2) 

NN=1 

540  DO  560  1=1.NC0N 
XS I )=0.0 
DO  550  J=1.NC0M 
550  X( I )=XI I )-B( I.J>*ZIJ) 

XI  I  )=X( I )+G( I ) 

560  X( I )=EXPF(X(  I ) ) 

564  IFIL)  565.565.610 


TABLE  7  -  Continued 


bbb  00  bib  J=UNC0M 
(jG(J)=0.0 
DO  b70  I  =  lfNCO.N 
5  70  GG(  J)  =G(i(  J)-i-bC(  I  ♦  J)*X(I  ) 
bib  GG( J)=GG( J)-G( J) 

DO  580  J=1*NC0M 
DO  580  K=i*NCOM 
A( JfK)=0,0 
DO  580  I  =  UNCON 

A(  J*K  )=A(  J*K)+dC(  I  ♦  J)*b(  i  tiCJ  I  ) 

^^80  AEQ(  JfX)=A(  J,K) 

ARG6=0 

LLL=I5IME0(  iO*NCO;'^,  1  ^AEQ  tCGf  ARG6  ♦Dh ) 
IF(LLL-l)  585*585f583 

583  PRINT  584 

584  F0RMAT(8H  STOP  5) 

585  ZZ=0.0 

DO  590  J=1*NC0M 
Z(J)=Z(J)+AEU(J*1) 

590  ZZ  =  ZZ+ABSF(AEO( J*l)  ) 

IF  (ZZ-CDZ)  600*600,540 
600  L=1 

GO  TO  540 
610  XN=0, 

DO  6Z0  J=1*NC0M 
UNblJ) 

620  WMM(J)=X(;) 

DO  630  I=1*NC0N 
630  XN*XN+BR( I)*X( D 
XN=1./XN 
WM=WMO/XN 
H2=0.0 

DO  640  I=1*NC0N 
640  H2  =  H2*»-H(  I  )*X(  I  ) 

H2=H2/WM 
705  CPRrO, 

DO  710  I=1*NC0N 
710  CPR=CPR+CP( I)*X( I) 

DO  730  I=1*NC0N 
DH( I )=0, 

DO  720  J=1*NC0M 
K=NB( JI 

720  DH( I )=DH( I )-6( I ♦J)»H(K) 

730  DH(  I  )=DH(  I  I  ) 

DO  750  J=1*NC0M 
GG( J)=0, 

DO  750  I=1*NC0N 

7  50  GG(  J)=GG(  J)+BC(  I  ♦J)<^X(I  )»0H(  I  )/(R<fT2) 
DO  752  I=1*NC0M 
DO  752  J=1*NC0M 
752  AEQ( I  ♦J)=A( I fJ) 

ARG6=0 

LLL=ISIMFO( 10*NCOM* l*AEQ*GG*ARGb*DZ) 
IF(LLI-l)  760*760*758 

758  PRINT  759 

759  F0RMAT(7H  STOP  5) 

CALL  EXIT 

760  DO  775  I=1*NC0N 
BETA(  I)  =  0. 

DO  774  J=1*NC0M 


TABLE  7  -  Contined 


7  74  RrTA(  I )=RFTA( I )-R{ I ,J)*ALQ( Jtl  ) 

11b  RF  TA(  I  )=PFTA(  D+DHC  ])/{R*12) 

CPR=CPR/WM 

CPE=CPR 

PLNMsO* 

DO  777  IslfNCON 

7  77  PLNM=PLNM-WMP( I ) *X ( I ) *bt  F A ( 1  ) 

PLNMsPLNM/ ( T2*WM#*2 ) 

00  780  !=lfNCON 

780  CP£*CPE-»-H(  I  )*X(  I  )#8ETA(  1  )  /  ( WM#  T2  ) -fHC  n*X(  I  )^^PL^|M 
CP2=CPE 
RETURN 
END 


TABLE  8.  THERMODYNAMIC  DATA:  COEFFICIENTS 
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.5206998  -1.76641330  L 179612700  -.381344800,060163981  -.00367756  28790,13900  -.59613732 


/ 

TABLE  9.  Sample  Input  Data  for  Computer  Program  for  Mixture  Initially 

at  60^F  and  760  mmHg* 


Data  for  Spec(a)  Subroutine 


1.  288.72  .000001  .OOOOJ  20.  lt». 

10. 

6  2  2 
12 
3  2 

1.7747932  12.496609  -33.09369940.68232/  -23.2340093.038 /d62  -b63.d.00683# Oi  >  /  i /ci 
3.3963743  -.304322673.4248489  -3.4348 /62i‘.233V84b  i  1 996243- 1023.48 7^:j. 2  i  i 

4.0821991  -1.11136904ji0978V67  -3.0381%0li«i289Hib  -.  lb:i28Hol-302(J^o0i->-  . JdoooO  >2 
4.0067420  -2.41231144.3443704  -4.07061442.0392313  -•44010 ? 3333/4.8 12 9  -.i/brv./i 
2*4999243  23472*/bl  — .43:y4o3‘^3 

3.0439461  -2.25983073.9647901  -3. 4822i0ul. 304223  /  -.2341 /31  / >91 34.002  2.33d3v4^ 
2.8699448  .  /2232810  .06 3 994338-. 09u  72232 .02 J4H68b3 -.00 14  /b82  '  /4l . 1 74  /x-. 3  /Od / 43 j 
3.8411221  -3.04038982. 248694:>  -• /2i662G8.ilUl8ulC  -.•JOb32d.:b-2228.27b2-d 7 o 
8.U262346  -6.4726129D.3982138  -1 .9863403.32330238  -.01997833-3229  /.  791-23.  34.'*2 
4.S902840  -1.71685631.6660167  -.60*76663  7 . Iu002323  x988  #3223.8  lb  /  -2.bb4H^*d/’ 

2.4999243  23472./61  -•h9>40:^:#3 

3.5206998  -1. 76641331.1796127  381 34480. 06016398 1-. 0036 7  7 362d790. 139  -.39bi3  /  i>2 

18.016  32. 

1.  -.5 

1. 

1. 

.3  .25 

.5  -.25 

.5 

^2.  1. 

H2  02  H20  OH  H  0 

.45  .08 


Data  for  Spec(b)  Subroutine 


1.  288272  .000001  .OJOOl  Xo.  lO. 

2  2  2 
12 
12 

1.7747952  12.496609  -33.09569940.68232  /  -23.2:>40093.03d7od2  -oo3.dG086  3 •  0.1  3  V 7>o 
3.3965745  -.304322673.4248489  -3.4348 762 1 .2 359846  -. 1 x996243-i023.4b725.2 Ti 3^ >o 
2.8699448  .72232810  . 063994338-. 09072252 .020448865-. 0014  7682- /^l .  1  >4 /!-.:>  706 
5.8411221  -3.04058982.2486945  -.72166208.11018010  -.00652828-2228.2782-8.9 7b5j>d5 
2.016  32. 

1. 

1. 

2.  1. 

H2  02 

.666 

5. 


.333 


TABLE  9  -  Continued 


Data  for  Spec(c)  Subroutine 


!• 


288.72 


.000001 


rOCOOl 


1. 


10, 


10« 


4C 

1 

3 


6 

2 

2 


1.7747952  12.4966^9  -33.095690 
3.3965745  -. 304322673.4248489 
4«Ci821991  -1.11136904.0978967 
4.CC67420  -2.41251144.5443704 
2.4999245 

3.0439461  -2.25985073,9647901 
2.0699448  .72232010  .">63994338 
5.8411221  -3.04058902.2486945 
8.0282346  -6.47261295.5982138 
4.3902840  -1.71685631.6660167 
2.4999245 

3.5206998  -1.76641331 


40.687327  -23.2540^95. 
-3.43487621#2369846  -. 
-3.05815011.1289416  -. 
-4.07C81442.C592515  -. 

-3.48223001.5042257  -. 
-.090-»2252. 020446066-. 
-272166208.11018010  -• 
-1.9865405.32350238  -. 
-.60766657.10002323  -. 


.1796127  -.36134480.06016398]-. 


0387882  -863.800883 
11996245-1023.48725 
18128481-30283.013- 
440107533574.8129  - 
25472.761  - 
2541731729134.002  2 
00147682-741.17471- 
00652828-2228.2782- 
01995635-32297.793- 
006198873223.8187  - 
25472.761  - 
0036775628790.139  - 


.0137796 

.2713998 

•38880032 

•17879530 

.45946595 

.5385942 

•57067440 

8.9783483 

23.324212 

2.8644487 

•45946595 

•59613732 


18.016 

!• 

1. 

•  5 
.5 


2. 


.45 


h; 


32. 

-.5 

1. 

.25 

-.25 

.5 

1. 

.C8 


02 


H20 


OH 


i 

I 

} 

] 


NOMENCLATURE  FOR  TABLES  10.  II.  and  12 


Initial  Condition  of  Mixture 

PI  ~  ‘^1 ,  o.'tm 
TEMP  -  T  ,  "K 
ENTHALPY -y, 

SD  VEL  -  a.,,  mA«c 
VI  -  V.,  cmVgm 
CPI  -  Cf»i/R. 

GAMMA  1  -  Vi 
WMl  -  m 
ENTROPY  -  S./R, 

Properties  Downstream  of  Wave 

PRESS  R  -  P 
TEMPR  -  © 

VEL  R  -  V 
MACH  1  -  Mi 
EMACH  2  -  M,,2 
F  MACH  2  '  M(  2 
ESDVEL  2  — 

FSDVEL  2  ,  m/sec 

F  GAMM  2  - 
E  GAMM  2  -  >< 

FSPHT  - 
ESPHT  -  C^Ja, 

ENTH  2  -  ^ 

MOLWT  -  mz 
ENTROPY-  s//?, 


Table  10: 


Sample  Output: 


Hugonlot  Curve  Computed  Using  Spec(a) 
Subroutine  and  Input  Data  In  Table  9 


HUCpMICf  ft«  USrCUS  PllTb«tS»*MO»U.tlOM  CIrliitlilCS  WklVMSlI* 

INIIIH  CCMClTICk  C*  riltult 


pi»  e.iccscf  .itu 

04C*»/Cf 

(«TSGi>y«  tV.IOII 


ZU.IIjO  CESK 
Cl»l«  J.*I4 


(•dIhUKV*  '<i.U 

uwn*  i.4e4 


:c  Ml*  it*.  u*u>/%u 

kfl*  I2.CI1 


cokoosiTici  cr  kuTiie 

M2  •  l.ttttTt  CJ 
a2  •  £.)I)IIE*C3 


tMiSHcctkirtc  mriatits  cn  tmi  Hucciiot  ciRvf 


OOfSS  1 

ler*  * 

871  *  MCMl 

M4C7>2  (M4L62 

f$L'7(l2 

7109(12 

8C4882 

(6*882 

89887 

«8m7 

11762 

6U67 

7..70C87 

10.0300 

12.1112 

0.9*4)12.902711. 319911. 4 74( 

11S4.7  * 

1444. *2) 

1.21*7 

1.212) 

4.99)) 

24.18)9 

*.699 

14.*;) 

*9.2279 

h 

11.0300 

12.22*9 

C,4l:9  1.9194 

2.79*4 

2.8064 

1961.. ^99 

liC3.7J) 

1.21*8 

1.2114 

4.99*7 

24.0914 

4.4)9 

)..64& 

29.19.3 

12.33C3 

12.1321 

0.*4u7  7.312* 

1.94}t 

2.971) 

1966.474 

I9S6.997 

1.2149 

1.2149 

4.*;t3 

2). 9449 

1C..M 

14.646 

2. .17- 7 

.1 

ll.OOCb 

12.172C: 

0.7*19  6.2916 

1.44*9 

1.7123 

i:72.:7* 

>912.297 

1.217r 

1.2159 

4.*<)) 

2).097) 

1..971 

14,627 

I-*.:*.:* 

14.33CC 

12.417* 

C.71.6  9.8621 

1.41*2 

1.4447 

1977.911 

1917.74* 

1.217) 

1.21(0 

4. *3*. 

2).70** 

11.1)1 

14.5v6 

29.1)7i 

19.CM6 

12.9324 

0.4*64  9.6181 

1.2912 

1.144) 

19*2.861 

1921.^88 

1.2179 

1.2177 

4. *389 

2).  7290 

11.487 

14.>*1 

* 

14.03CC 

12.9640 

0.4474  9.9.98 

1.1894 

1.2199 

19*8. :«9 

1928. )j4 

1.2177 

1.2187 

4.4)16 

2).**93 

12.2)1 

14.9..4 

2..U  re 

17.03(3 

12.4214 

0.4111  9.428S 

1.1071 

1.19C2 
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COMPOSITION  BEHIND  SHOCK  WAVE 


TABLfE  13.  Computed  Properties  of  CJ  Detonations  in  Hydrogen- Oxygen  Mixtures. 

Initial  Cooposition:  ^2*^2 
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Fig.  1:'  Determination  of  a  Point  on  the  Hugoniot  Curve  in 

the  V-P  Plane 
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Fig.  2:  Determination  of  a  Point  on  the  Hugon.:.i 
theV-Q  Plane 
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Fig*  3‘  Diagram  Showing  Hugoniot  Curve i  Ranklne 

Hugonlot  Curve,  and  CJ  Rayllegh  Line  for 
2H2+O2  Mixture  at  60®F,  760  mm  Hg. 


IP  JWP'W  J.. 


Tt  CO  CM  T-  P 

ON-  HDVN 

% 


Fig.  4:  Determination  o£  the  CJ  State  in  tl^e^VPPl»“® 


GO  TO  SPEC.  WITH 
INITIAL  DATA. 


GO  TO  MAIN  WITH 
p  SPECIFIED 


SUBROUTINE  SPEC 

(TABLE  3.4.5  ) 

SPEOFY  VALUES  OF  p  AT 
WHICH  HUGONIOT  CURVE  IS  TO 
BE  SOLVED  OR  FIND  B,  .ORf^j. 
PRINT  OUTPUT.  ^ 


SUBROUTINE  MAIN 

(TABLES) 


EVALUATES  THE  THERMODYlvlAMIC 
PROPERTIES  ON  THE  HUGONOT 
CURVE  AT  THE  QVEN  p. 


GO  TO  BPiNKWlTH 
p  AND  T  SPBCIRED 


SUBROUTINE  BRINK 

(TABLE?) 


EVALUATES''EOUILIBRIUM  AT 
SPECIRED  p  AND  T. 


RETURN  FOR  MORE 
DATA  OR  EXIT  FROM 
PROGRAM. 


RETURN  TO  SPEC 
WITH  SOLUTION  TO 
HUGONIOT  CURVE 
AT  SPEDFED  p. 


RETURN  TO  MAIN 
WITH  EOUIUBRIUM 
COMPOSITION  AT 
GIVEN  p  AND  T. 
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Fig.  7:  Hiigoniot  Curves  for  Hydrogen-Oxygen  Mixtures 
Initially  at  60^F 
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Fig.  9:  Hugoniot  Curvet  for  Hydrogen- Oxygen  Mixtures 
InitUlly  at  -180®F 


Fig.  10:  Hugoniot  Curve  in  Dimensionless  Temperature -Entropy  Diag 
for Mixture  Initially  at  60^F,  760  mmHg. 
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Fig.  12:  Influence  of  Initial  Pressure  on  CJ  Pressure  Ratio 
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Fig.  13:  Influence  of  Initial  Temperature  on  CJ  Pressure  Ratio 
forBH^'^C^Mixture 
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Fig.  14:  Influence  of  Initial  Pressure  on  CJ  Temperature  Ratio 
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Fig.  15:  Influence  of  Initial  Pressure  on  CJ  Detonation  Velocity 
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Fig.  16:  Influence  of  Initial  Temperature  on  CJ  Detonation 
Velocity  for  +  Mixture 
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Fig,  18:  Influence  of  Initial  Pressure  on  Equilibrium  Composition  of 
Product  Gases  for  3W^■^Op  ^  Iritially  at  ZOO^F 
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Fig.  20:  Influence  of  Initial  Pressure  on  Equilibrium  Composition  of 
Product  Gases  for  Mixture  Initially  at  -50^F 
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Fig*  22:  Influence  of  Initial  Pressure  on  Equilibrium  Composition  of 
Product  Gases  for  214^ Mixture  Initially  at  60°F 
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Fig.  24:  Influence  of  Initial  Pressure  on  Equilibrium  Composition  of 
Product  Gases  for  Mixture  Initially  at  60®F 
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■Fig.  25: 


Influence  of  Initial  Pressure  on 

Product  Gases  for 


Equilibrium  Composition  pf 
InitL?lly  at  -180  F 
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Fig.  26:  Influence  of  Initial  Pressure  on  Specific  Heat  at  the  Constant 
Pressure  of  Products  at  CJ  State.  Initial  Temperature  60®F 
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Fig.  27:  Influence  of  Initial  Pressure  on  Specific  Heat  at  the  Constant 

Pressure  of  Products  at  CJ  State.  Initial  Temperature  -ISO^F 
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Fig.  31:  Influence  of  Initial  Pressure  on  Von  Neumann  Spike  Velocity  Ratio 


